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ABSTRACT 


This thesis examines the performance of SRI's Ambient Communications 
(A.MBCO.M) model for high-latitude propagation prediction. It is one in a series of 
studies, conducted at the Naval Postgraduate School, to establish the relative merits of 
several computer-based propagation prediction models using a standard set of measured 
data. 

AMBCOM modeled the propagation path between a transmitter located in the polar 
cap region and several midlatitude receiver sites. Model predictions were matched to 
measured data obtained during two high-latitude communication experiments (cam¬ 
paigns). The absolute dilTercnce between model signal-to-noise ratio (SNR) and meas¬ 
ured SNR was considered as error. Error statistics were accumulated to show the 
distribution of the error by campaign and frequency. The percentage, by frequency, of 
matched AMBCO.VI predictions in reference to total predictions for a given frequency 
was considered a measure of AMBCOM performance. 

A.MBCO.M exhibited small absolute values of average error, i.e., 7-11 dB, and high 
percentages of matched records. The average error was typically distributed between -20 
and +20 dB. Unfortunately, these are only relative measures of model performance. 
The site antenna and environmental data used to model high latitude campaigns were 
estimated not measured, and some variation in AMBCOM results may be attributable 
to poor estimates. The measured data were not designed specifically for model vali¬ 
dation, and further comparisons are needed with new measured data. 
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I. INTRODUCTION 


A. GENERAL 

Duriiig the summer of 1986 L'RSl sponsored a Workshop on the Practical Aspects 
of Ionospheric Modeling. The conference was attended by thirty-three scientists from 
five countries. The goal of the conference was to bring together ionospheric scientists 
and practioners to discuss ionospheric modeling. One of the recommendations which 
resulted from this conference was the following: 

There is a need to have a standardized set of data that can he used to verify 
ionospheric models. These data should be independent of the models that are being 
verified. The data need to he representative of the entire range of season, solar, and 
geomagnetic activity variations that is observed in practice. [Ref 1; p. 7-13] 

The purpose of this thesis is to examine the accuracy of propagation predictions 
from the Ambient Communications (.-W'''M i computer-based model. This is done 
by comparing the signa!-to-noisc ratio (S’...) data generated by .AMBCOM with the 
SNR data recorded during Project N()N( LNTRK.'. 1 he Project NONCENTRIC data 
are the signal and noise strengths of transmissions which were broadcast from one site 
via an omni-directional antenna to several receiver sites during three high latitude com¬ 
munication experiments (campaigns) [Ref 2). 

This thesis is divided into five chapters. The first chapter is a brief overview of the 
morphology of the high latitude ionosphere as it applies to A.MBCO.M. Chapter one 
also describes previous computer-based model evaluation work. The second chapter will 
describe the A.MBCOM algorithms. Chapter three will characterize the Project 
NONCENTRIC data and describe how the data were assembled. The fourth chapter 
will compare the AMBCOM data with the Project NONCENTRIC data. The thesis 
then concludes with recommendations in chapter five. 

B. THE POLAR IONOSPHERE 

1. The Earth's Magnetic Field 

The earth's magnetic field may be approximated by an earth-centered dipole 
field which is tilted slightly out of the axis of the earth's rotation. At the present time 
the northern pole of the dipole is located approximately at 78.8° N, 70.9°W using the 
geographic coordinate system [Ref. 3: p. 83]. As illustrated in Figure 1 on page 2 the 
earth's magnetic field is distorted from the expected dipole shape by the effect of the 


1 



solar wind. Traveling at speeds near the earth of 200 to 700 kin/s, the solar wind bends 
the earth's magnetic field into a teardrop shape with the tail of the field located opposite 
of the earth's sunward side [Ref. 1: pp. 2-ld;2-21]. I he tail of the field can extend out¬ 
ward several earth radii depending upon the force of the solar wind [Ref !: p. 2-17], 



Figure 1. The cartli'.s magnetic field shoniiig the effects of the solar wind. [From 
Ref. 6: p, 2-25] 


AiMBCOM uses the three-hour A', index of worldwide magnetic disturbance to 
specify the current state of the earth's magnetic field, fliis index is based upon local K 
indices which are semi-logarithmic values prepared at twelve selected observatories lo¬ 
cated between 48“ and 63“ geomagnetic latitude, north and south. The local K indices 
are values assigned by a local observatory to describe the condition of the planetary 
magnetic field at the measuring site. 1 hese local K values are then corrected and used 
to calculate the planetary A, index. (Ref 4; pp. 49-51] 

The A', index is calculated at three hour intervals for a total of eight periods per 
day. The A, index ranges in value from zero, the least disturbed state, to nine which 
represents the most disturbed magnetic field, lypically the A, indices are calculated 
monthly and published in tabular form. The A, indices used for this thesis were pre¬ 
pared by the Geophysikalisches Institut, University of Gottingen and tabulated in the 
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Journal of Geophy sical Research (Ref. 3]. These indices were as low as zero and as 
high as seven which means that the geomagnetic field varied from a nondisturbed to a 
considerably disturbed state. 

When studying or modeling the ionosphere, the geomagnetic coordinate system 
is commonly used to map the ionosphere to an earth-bound coordinate system. This is 
understandable in light of the fact that the earth s magnetic field plays a prominent part 
in the formation and variability of the ionosphere. The geomagnetic coordinate system 
is based upon the geographic location of the earth's magnetic poles. The geomagnetic 
equator is displaced 90° from the geomagnetic poles. The longitudinal origin for this 
geomagnetic system is the meridian line which passes through the north and south 
geomagnetic poles and through the geographic south pole [Ref 3: p. 84]. Davies (Ref 
4; p. 40] relates the two systems with the equations 

sin O = sin (/) sin 0o + cos 0 cos 0o cos(/ — /q) (I) 


and 


sin A = 


( cos 0 sinCz — >.(,)) 
cos <1> 


( 2 ) 


where 

00 = geographical latitude for the northern geomagnetic pole, 

/o = geographical longitude for the northern geomagnetic pole, 

0 = geographical latitude, 

). = geographical longitude, 

O = geomagnetic latitude and 
A = geomagnetic longitude. 

Although A.MBCOM makes the geographic-to-geomagnetic coordinate conversion for 
the user, an understanding of the geomagnetic coordinate system is useful with regard 
to understanding the AMBCOM model. 

2. Particle Precipitation 

As the solar wind passes the earth, electrons, protons and other particles are 
captured by the earth's magnetic field. A captured particle does not fall directly to the 
e:’ rth but follows a spiral path along the magnetic field lines with velocity components 
parallel and perpendicular to the magnetic field. As a particle moves to higher latitudes 
the magnitudes of the components of its velocity vector change, but the total particle 
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energy remains constant. This is due to the bending of the magnetic field lines relative 
to a particle's velocity as the earth's magnetic field returns to the ground at the magnetic 
poles. Increasingly the particle's velocity becomes perpendicular to the magnetic field. 
The forward motion of the particle towards the pole stops and then reverses. The par¬ 
ticle then begins follow a spiral path towards the opposite geomagnetic pole, and, again, 
the components of the particle's velocity vector change as the magnetic field lines bend 
earthward. Forward motion again stops and reverses. This point of reflection is called 
the mirror point. 

Depending upon the location of the mirror point, the particle can be trapped in 
the earth's magnetic field indefinitely. If the particle's mirror point, which is dependent 
upon its pitch angle when it entered the earth's geomagnetic field, is above the atmos¬ 
phere then the particle will continue to reflect between the north and south mirror 
points. If the mirror point is within the atmosphere then the particle may collide with 
other particles in the upper atmosphere, lose energy and precipitate into the polar 
ionosphere. This particle precipitation leads to increased radio signal absorption, i.e. 
auroral absorption and polar cap absorption (PCA), as well as visual displays of aurora 
borealis. Energetic electrons are the major cause of auroral absorption and the visible 
aurora. These particles have an energy level of about 10-100 KeV. PCA is attributed 
to the heavier but more energetic 10 MeV protons. [Ref. 6: pp. 2-18;2-20] 

3. The Auroral Oval 

The auroral oval, see Figure 2 on page 6, is the region where the visible aurora 
is most often observed (between 64° and 70° north geomagnetic latitude) [Ref 4: p. 51]. 
The visible aurora is the result of the precipitation into the upper atmosphere of 
electrons which have been carried to the polar regions by the earth's magnetic field. 
Since the amount of precipitation depends upon the number of particles available, any 
occurrence which increases the number of particles captured by the earth's magnetic field 
will increase particle precipitation in the auroral zone. During a solar event, such as a 
solar flare, large numbers of high energy particles are ejected from the sun. These par¬ 
ticles speed towards the earth and cause a disturbance in the earth's geomagnetic field 
which leads to particle precipitation in the polar regions. Jones [Ref 2] makes the 
metaphorical comparison of particles traveling along the magnetic field lines as being in 
tubes something like a tube of toothpaste. When a disturbance occurs in the earth s 
magnetic field the particles (which are moving between the mirror points) are "squeezed 
out" into the high latitude upper atmosphere. This results in higher propagating fre¬ 
quencies and auroral absorption. 
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The size and location of the auroral oval at any given time is variable because 
of its dependence upon particle precipitation which varies with solar activity . The oval 
is usually thicker and extends equatorward on its midnight side and migrates poleward 
on the noon side (Ref 7: p. 348]. The midnight side of the oval is where substorm ac¬ 
tivity begins before it follows the electrical potential gradient toward the oval's noon side 
[Ref 8: pp. 431-437], In addition to the visible aurora these particles also produce 
ionization which results in auroral absorption, sporadic E, and other high latitude phe¬ 
nomena which do not always coincide with the location of the auroral oval [Ref 4: pp. 
332-336]. 

Equatorward of the auroral oval is a region of low ionization density called the 
midlatitude trough. The midlatitude trough is caused by ionospheric currents resulting 
from the large electrical potential difference between the day and night side of the polar 
region. Hargreaves [Ref 1: p.7-7] describes the trough as a nighttime occurrence be¬ 
tween 60° and 65° geomagnetic latitude during all seasons. During the summer months 
the midlatitude trough occurs closer to midnight. At night and during magnetic storms 
the trough moves equatorward. The midlatitude trough is of interest to the commu¬ 
nicator because it effects propagation paths. Critical frequencies within the trough can 
be as low as 1 .MHz. Higher frequency transmissions along paths with reflection points 
in the midlatitude trough are not reflected downward and pass into space. 

4. High Latitude Propagation Anomalies 

The higher latitude ionosphere is not as dependent upon direct solar x-ray and 
ultraviolet radiation for the creation of ionization layers such as the D, E, sporadic L 
or F layers. The layering of the polar ionosphere is produced by a combination of solar 
radiation and particle precipitation. The contribution of particle precipitation is ob¬ 
servable as increased ionization density in the E layer and a reduction in f, ionization 
density during periods when the solar zenith angle (x) is greater than 90° [Ref 9: p. 28]. 
As can be seen in Figure 3 on page 7, auroral absorption exhibits two maximums, one 
at midnight and one at noon. Hargreaves [Ref 10: p. 1358] states that there is a definite 
correlation between the amount of magnetospheric disturbance and increased auroral 
absorption. An increase in the A", index wmuld presage an increase in auroral absorption 
due to particle precipitation into the auroral oval. Hargreaves [Ref 10 : p. 1358] states 
that increases in the density and speed of the solar wind correlate with changes in A”,. 
An increase in solar activity results in an increase in the number of particles captured 
by the earth's geomagnetic field. One result of this process is more auroral absorption. 
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Figure 2. The auroral oval and the propagation path from Cl3'de River to site O. 

In addition to tlie auroral absoiption, communication signals at high latitudes 
may (bllow non-gicat-circle (NCiC) propagation paths. At lower latitudes IIF signals 
generally follow a great circle path from transmitter to receiver. High latitude signals 
can be reflected in a direction that is perpendicular to the great circle propagation path, 
llunsuckcr (Ref. 7: pp. 348-351] states that NGC propagation may be attributable to 
sidcscattcring due to irregularities in the ionosphere. The incidence of NGC modes have 
been found to vary with the amount of magnetic activity. The final result for the com¬ 
munications engineer and the AMHCOM user is that signal strengths predicted by a 
model may not match measured signal strengths because the effects of NGC are not in¬ 
cluded in the model. During times when magnetic activity is liigh, the inability of the 
computer-based model to accurately predict signal strengths at given frequencies may 
be attributable in part to non-great-circle paths. 

I he occurrence of sporadic II {E,) at higher latitudes is important to communi¬ 
cation in the IIF band. E, provides additional paths for propagation at IIF frequencies 
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Figure 3. A typical day in auroral absorption. [From Ref. 10: p. 1361) 


as well as some attenuation of signals propagated by other modes. Ilunsucker [Ref. 7: 
p. 353] described the effect of E, resulting from measurements taken for propagation 
paths between fhule (Greenland) and College (Alaska) as well as for a path from 
Andoya (Norway) and College as follows: 

In particular, for winter-night sunspot minimum conditions on the Thule-Collcge 
and Andoya (Norway)-College path, E, modes had an average occurrence over 50%. 

E, MOF's as high as 46 MHz with typical values of 18 MIIz were observed. 

As can be seen in Figure 4 on page 8 the percentage of sporadic E varies with time and 
latitude. In the auroral oval the average percentage of E, is at a maximum during the 
evening hours and at its lowest at noon. Ilunsucker (Ref. 7: p. 351] states that E, is more 
prevalent and with higher maximum observable frequencies (MOF's) during the winter 
months compared to summer months. E, w'as included in some path calculations made 
by AMBCOM. In addition some effort was made to adapt the percentage of E, to the 
time of day and year. The results from E, calculations were then compared to data from 
AMBCOM which included no E, calculations. (Ref. 7: p. 351] 

The polar ionosphere does include F layers similar to those found at the 
midlatitudes. The distinction between the polar F2 layer and the FI layer may not be 
as dramatic as at low'cr latitudes. IT mode propagation is frequently observed during 
the summer, especially in the morning. 1 his is because the F2 layer critical frequency 
decreases with increasing latitude at a faster rate than the FI critical frequency. [Ref. 7: 
p. 335]. The FI mode is the dominant means of long-range propagation during periods 
of low sunspot activity, particularly during the low' end of the eleven year sunspot cycle. 


7 








Figure *4, Tlie average percentage of sporadic E by geographic latitude and time of 
day. (From Ref. 7: p. 352J 


Vondrak (Ref. 9; p. 38) noted that tlic ionization density profiles recorded by the 
Cliatanika radar showed a valley between the H layer and the F2 layer when the layer 
was not sunlit. This valley disappeared duiing the day and during periods of high particle 
precipitation. 

C. PREVIOUS PROPAGATION MODEL STUDIES 

Ilow accurately can high latitude propagation be predicted using the computer 
models that arc presently available? Presently there is not a complete answer to this 
question. Although individual studies have been carried out comparing the output from 
computer models to the data recorded from transpolar and polar communication paths, 
these studies were not intended to be comprehensive comparisons of the accuracy of 
numerous models relative to one another using one set of recorded transmissions. 
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1. The Rockwell International Study 

During fiscal year (FY) 1987, Rockwell International completed a transauroral 
high-frequency (HF) radio transmission experiment (TAHFE). The experiment con¬ 
sisted of recording, at Cedar Rapids, lA, transmissions from Barrow, AK. The trans¬ 
missions were made by stepping through a selected set of frequencies at specific times. 
Figure 5 on page 10 is a block diagram of the equipment used for this experiment. 

In addition to the transmissions, data were recorded which described the 
ionosphere, solar activity and HF channel parameters at any given time. 
Solar geophysical data were collected from the U.S. NOAA Space Environment Services 
Center (SESC) [Ref 11: p. 2]. The HF channel parameters where measured by Rockwell 
International engineers using the advanced link quality analyzer (ALQA). These data 
were recorded on computer disk for later analysis. Using SESC data and with the aid 
of solar event alerts from the Geophysical Institute at Fairbanks. AK, the experiments 
were able to match measured channel parameters to solar activity. [Ref 11: pp. 1-3] 

Three computer models, i.e., lOXCAP, .MlNl.MLT and A.MBCO.M, were used 
to generate predicted maximum useable frequencies (.VIUF) using the recorded solar- 
geophysical data as a primary input for the prediction models. All models produced 
good approximations of MUF values that were actually recorded during the experiment. 
AMBCO.Vl was able to predict the effects of solar events upon the transpolar path such 
as the a communications blackout on 12 November 1986. [Ref 11] 

2. Hiah-Latifude HF Predictions From RADAR C and AMBCOM 

Nikhil Dave of the Naval Ocean Systems Center (NOSC), San Diego, 
California, conducted a comparison of the data produced by two models (R.ADAR C 
and A.MBCOM) and polar sounder data recorded for the Winipeg-Resolute Bay path 
during 1959 and data recorded during 1964 over the two paths (Andoya-Ft. Monmouth 
and Andoya-College during 1964). The former data were recorded during a period of 
high sunspot activity and the latter data during a period of low activity. The intent of 
the study was to produce a predicted MUF and to compare it with recorded data. 
AMBCO.M was found to be more accurate than RADAR C because AMBCO.M in¬ 
cludes adjustments for auroral absorption and other polar anomalies. [Ref 12 ] 

3. The Naval Postgraduate School Studies 

The work at the Naval Postgraduate School was intended to compare the per¬ 
formance of as manv models as possible to the same measured data, antenna gain, and 
control parameters. NPS students exantined the relative performance of lONCAP, 
PROPHET, and A.MBCO.M. In the future, a new high-latitude propagation model, 
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Unmanned Remote Site ^ 
Barrow AK 


4750 KM 


Control Site 
Cedar Rapids lA 


l iguie 5. TAHFE cquipiiiciit suite tiiat was used for tlie Rockwell experiments. 
|I roni Ref. 11] 


ICEPAC, will be studied. Ihe comparison of data generated by lONCAP and 
PROPHET were completed and published [Ref. 13 and 14 ], and may be compared to 
this A.MBCOM analysis (see Appendix D). It is hoped that the results of these studies 
will guide the future creation of a standard model validation data base and the estab¬ 
lishment of guidelines for computer-based propagation prediction model evaluation. 
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II. THE AMBIENT COMMUNICATIONS MODEL (AMBCOM) 

A. INTRODUCTION 

This chapter provides a brief description of applicable portions of AMBCOM. De¬ 
tails of the internal functioning of AMBCOM are contained in the user's guides pub¬ 
lished by SRI International [Ref. 15,16 and 17j. 

AMBCOM is designed for batch processing using card images as input. Separate 
programs support each of the function areas modelled, i.e., the ionosphere model is 
generated by a program which passes its data to another program for raytracing calcu¬ 
lations, The system flow is first explained followed by the programs and inputs which 
generate the high latitude data. 

B. SYSTEM FLOW 

The AMBCO.M system is a multiprogram batch system written in FORTR,A\ and 
was run on a VAX 3100 workstation at NTS. Programs NATGEN, RA’lTR/X, and 
COMEFF were used to generate the high-latitude data, as shown in Figure 6 on page 
12. Data are passed between these programs by means of saved data files. [Ref 15 : 
pp, 107-1551 

Since A.MBCO.M programs pass data through standard files, the user must be care¬ 
ful with regard to the timing of program execution. Several campaigns were modeled in 
this study. If two campaign decks were executed simultaneously and the programs did 
not terminate due to file access collisions, then the accuracy of the contents of a passed 
file could not be guaranteed. The user has the option of either changing the names of 
the passed files so that these file names differ from one campaign input stream to an¬ 
other or to simply run one campaign input stream at a time. The latter method was 
chosen so wiiat each execution of AMBCO.M would not share the assets of the VAX 
3100. 

The execution input stream were constructed for all 25 days of a particular cam¬ 
paign in only one computer run. Separate AMBCO.M input streams were created for 
each site during the campaign. This entails the construction of an input stream for 
NATGEN, R.^YTRA, and COMEFF, representing each three-hour interval and 
eight intervals per day for 25 days. Each campaign input stream was approximately 
13,300 lines in length, and required two and half hours of execution time. 


II 







Figure 6. AMBCOM system flow showing the passage of data between programs. 


1 he data produced by AMBCOM for each of the three hour intervals was accumu¬ 
lated’in files for later processing. This is not a problem because the VAX 3100 work¬ 
station allows the use of multiple generation data sets. These data sets were easily 
concatenated at the end of the campaign's e.xecution stream into one large file repres¬ 
enting all the data for one campaign. 

C. NA7GEN 
1. Overview 

The purpose of NATO EN is to model the ionosphere along the commudication 
path between two points. NATGEN builds a model of the F2, FI, and E layers at 
control points along the path. I he control points are evenly spaced at 100 km incre¬ 
ments with a maximum of 41 control points for paths longer than 4000 km. The control 
points are not necessarily located at the signal reflection points since raytracing is per¬ 
formed in RAYTRA. [Ref 15: pp.23-24J 

NATGEN begins by reading ionospheric coefficients provided by the Institute 
for Telecommunications Sciences (ITS) in Boulder, Colorado, based upon the month, 
day, and the current sunspot number. AMBCOM uses the ITS Blue Deckl as the 

I 

i A reference to earlier coefTicient files that were issued on color tiimmed compute!'cards. 
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starting point for the modeling of the Fj, F, , and E layer at each control point. Addi¬ 
tional parameters describing atmospheric noise, ground conductivity, and ground 
permitivity along the path are read and passed to RA^TRA. [Ref 15 : p.23] 

The semithicknesses, heights of maximum ionization, and the critical frequencies 
for the F^, F, and E layers at each of the control points are passed as outputs Irom 
NATGEN. ITS upper, median and lower decile values for sporadic E critical frequencies 
are passed to RAYTRA which performs all £, calculations. NATGEN also passes on 
the location of the transnutter and receiver sites, the time of year, time of day, the cur¬ 
rent index, the sunspot number, distance between control points, the number of 
control points, and the path length. (Ref 16: pp. 123-124] 

2. Layer modeling 

The ITS provides the vertical incident critical frequency (/„) for the E layer, the 
Fj layer, and sporadic E. The E layer critical frequency {f^E) is always set to the median 
decile values. For the sporadic E and the Fj 

layer, the upper, median, and lower decile values represent the critical frequencies 90%, 
50°o, and 10%, respectively, of the days for a given time and month [Ref 18: p.89]. 
N.ATGEN allows the user to choose which of the Fj layer values will be used as the Fj 
critical frequency (4F2) for all ionospheric calculations. The default f,F2 value is the 
median decile number. Two additional parameters for the Fj layer are the ratio of the 
semithickness of the layer (y„) to the height of maximum density (hj and the maximum 
usable frequency for a 3000 km path (M3000). [Ref 15: p.25] 

Using these coefficients, A.MBCOM models the ionosphere at the control 
points. The Fj, F, and E layers are represented as three parabolic layers with the values 
for the F, layer derived from the Fj and E layer models. The height and semithickness 
of the E layer are set at 115 km and 25 km respectively. The Fjlayer height is calculated 
using the ITS coefficients in a two step process. First the peak height of the Fj layer ( 
//F^j is calculated using the Shimazaki equation [Ref 15: p. 26] 


HPn = 


1490 

A/3000 


- 176. 


(3) 


This value is corrected for signal retardation with a height factor (Ah) equation [Ref 
15: p.26] 


Ah = 






y'm' 


(4) 
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where 


Ah = the height error, 

/ = the critical frequency of the layer, 
f = the transmitted frequency and 
y„ = the semithickness of the layer. 

The Ah factor is then subtracted from HP^ to produce a corrected height for the F 2 layer. 

The layer parameters are not represented by ITS data but are calculated using 
the E and layer parameters. The E, layer must overlap the layer by half of its own 
semithickness with the bottom of the E, layer set at 130 km. The critical frequency of 
the E, layer is calculated based upon the critical frequency of the E layer In the 

event that the F^ layer critical frequency is lower than that for E, 

/„E1 is set to 0.695 /„E2. [Ref 15: p.27] 

3. High-Latitude Corrections 

Before any correction for high-latitude propagation paths can be made 
N.-XTGEN must define the auroral oval using the Miller-Gibbs model [Ref 15: pp. 
28-30] with a value for the width of the oval in the path (X) and the location of the lower 
boundary of the oval (</>,,). The lower boundary is 

0^ = 71.9* - 2.5Kp - 5.1 cos[-|^ (t, - I)] (5) 

and must lie between the boundaries 

= 68.9* - a; - 5.1 cos[ ^ (r, - 1)] 

and 

0max = 70.9* -Ap - 5.1 cos[ -g- - 1)], (7) 

where 

t„ — geomagnetic time. 

The width is calculated as 

A'=7.0*-A^. (8) 
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The oval width is limited to between four and six degrees of latitude. 

XATGEN is designed to accommodate propagation modeling at high latitudes 
by correcting critical frequencies for the /y layer if„F2) and the E layer {f,E}. Corrections 
are made lof,Fl if the geographic latitude of the control point meets one of the follow¬ 
ing criteria; 

• The point is located at greater than 45° North (corrected geomagnetic latitude 

Uj)- 

• The point is located within the auroral oval. 

• The point is in the nighttime midlalitude trough. 

All the points used in this thesis are north of 45°N latitude and so the first criteria always 
applies. The correction applied to f^F2 at high latitudes is a two step process. First 
f„Fl is adjusted according to the first criteria above using 

foFl^=foFllTsi^-^Kp'\, ( 9 ) 


where 

f,f^A — adjusted for high latitude, 

= the F2 layer critical frequency from the ITS file and 
Sfr^ is a function of K^. 

All other corrections are applied to this adjusted value. [Ref. 15: p. 29] 

The second step is, to adjust/,F2 if the control point is in the midlatitude trough 
or the auroral oval. All control points that are north of the equatont'ard edge of the 
auroral oval are considered auroral oval corrections. Auroral oval adjustments are cal¬ 
culated based upon the location of the control point in reference to (f)^ and X. The 
correction is calculated using the Miller-Gibbs equation [Ref 15: p. 29] 

f,F2 = f,F2^ [l +0.4946 , (10) 


where 




( 11 ) 


This last parameter is the normalized distance referenced to the equatorw'ard edge of the 
oval. [Ref 15; pp.29-30] 


15 









The correction for the midlatitude trough accounts for the effect of this region 
onf^Fl. In this region the critical frequency can be as low as 1 to 2 MHz and a signal 
which enters this part of the ionosphere may penetrate into space. NATGEN makes 
corrections for the midlatitude trough only when the ionosphere is not sunlit, e.g., the 
solar zenith angle {y) is greater than 90° and the corrected geomagnetic time is between 
1800 and 0600. Reference 15, pages 30 thru 31, provides complete details of the cor¬ 
rection equations used by NATGEN. 

The high-latitude correction to f^E is much different than that for f^Fl. This is 
due to the fact that E layer propagation at higher latitudes is dependent upon solar ra¬ 
diation and particle precipitation. The location of this auroral E zone will vary with the 
time of day and with magnetic activity. NATGEN compensates for this by first model¬ 
ing the boundary of the auroral E zone and then calculating f E based upon the separate 
values for the E layer critical frequency due to 'ola. radiation (f^E(solar)) and particle 

precipitation (4£(awrora/))- 15; p. 31-35] 

The local K index is an indic.Uor of r>^rticle precipitation which results in higher 
E layer ionization densities. The magnetic field at higher latitudes is much different from 
that at midlatitudes and is not well represented by K^,. The index reporting stations 
are mostly located in midlatitudes and therefore the index is weighted more heavily in 
favor of midlatitude values. Hatfield [Ref. 15: p.32-33] states that does not produce 
high latitude modeling results which correlate weh with measured data, but using the 
local K index value does produce a stronger correlation to measured data. NATGEN 
has been adapted to estimate the local K index based upon the corrected geomagnetic 
latitude, the corrected geomagnetic time and the current index. This local K index 
value is used to calculate y^£(aurora/). 

The final step of the process is to combine the value iox f„E{solar) with that for 
f^Eiauroral). .According to Hatfield [Ref. 15: p.35] the final value fory^f is given by the 
equation 

f,E = ^^\.5f,E{aurora[f +m^olarf . (12) 

The f„E{aurorat) is given more weight than since particle precipitation heavily 

influ. ces the E layer at high latitudes. 
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4. NATGEN Input Variables 

AMBCOM is a batch system in which the input stream which controls the exe¬ 
cution of the system resembles a series of 80 column computer cards. A.MBCO.M was 
originally designed during the 1970s when card input systems still predominated [Ref 
15: p. 15]. Figure 7 on page 18is an example of statements used to execute NATGEN. 
The ASSIGN statements provide NATGEN access to the ITS file (ESSABLL'.DAT), 
the cnnversion tables for the geomagnetic coordinate system (RAGCOT.DAT), two 
output files (IONOS.DAT and NOISE.DAT), and execution control cards one through 
four. 

A.MBCO.M provides a large amount of flexibility in the control of data pro¬ 
duction. NATGEN models the ionosphere along any propagation path at any hourly 
time increment specified by the numbered control cards. The AMBCOM user defines 
the problem by specifying the transmitter and receiver locations, the sunspot number, 
the index, the time increment, and the use of auroral layers. 

Control parameters were chosen to model the ionosphere, as closely as possible, 
to the Project NONCENTRIC transmission paths. The intent was to execute the model 
from the perspective of a communicator who is attempting to estimate the possibility 
of communicating with another HF site under a given set of circumstances. Although 
the communicator does not have the advantage of using the sunspot number and in¬ 
dex tables for a current month, he should have a rough estimate of the current sunspot 
number and K indices at any given time. 

Figure 7 on page 18 show's the four control cards used to execute NATGEN. 
Card one is used to specify the use of auroral particle densities for all layer calculations, 
and card two causes NATGEN to run the ionspheric generator program. Card three 
defines the geographic latitude and longitude of the transmitter and receiver; negative 
numbers indicate a w'est longitude. Following the location information are the year, 
month, sunspot number, and A, index. Card four indicates the beginning hour, ending 
hour and time increment. 

Input streams w’ere divided into three hour periods over a 25 day period. This 
was for modeling the ionosphere during each of the three hour index periods 
throughout the day. The values for sunspot number and Aj, index were taken from the 
Journal of Geophysical Research [Ref 5, 19, 20, and 21]. The days of the month and 
hourly time increments were derived from Project NONCENTRIC data. 
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$DEL OWINTEKP.DAT;* 

$DEL ajMOUT.DAT;* 

$1345678901234567890123456789012345678901234567890123456789012345678901234567890 

$1 

$1 FREQUENCY MANAGEMENT BETWEEN A SPECIFIED TRANSMITTER AND RECEIVER 
$118 JANUARY 1989 

$l****** BEGIN KP PERIOD 1 •**•»*»♦**•*•«••*•«**»*•»•«***•*♦»*••****•••*******»• 

$1 *‘**INPUT SETUP FOR NATGEN**** 

$1 

$SET DEFAULT DKAO: [ AMBCOM | 

$ASSIGN DKAO: IAMBCOMIESSABLO.DAT FOROOl: 

$ASSIGN DKAO: IAMBC0M1RAGC0T.DAT FOR075: 

$ASSIGN DKAO:IAMBCOM)IONOS.DAT FOR002: 

$ASSIGN DKAO:(AMBCOMlNOISE.DAT FOR004: 

$ASS1GN SYS$INPlfr: FOR058: 

$ASS1GN DKAO: [WILSON)NATOOT.DAT FOR059: 

$R DKAO: (AMBCOM[NATGEN 

1 1 0 0 0 0 0 

BEGIN THESIS WORK WITH AMBCOM 

2 1 1 

3 70.00 -69.00 53.00 -01.00 1989 01.0 231.0 4.0 

4 0. 2. 1. 


y 


Figure 7. NAI GEN coiUiol cards u.sed to mode] h inter site O propagation data . 

D. RAYTRA 
1. Overview 

Tlie laytracing program (lOW'I RA) model is divided into two modes, point- 
to-point and radar. 1 he point-to-point mode performs raytracing from a transmitter to 
a receiver and was used for this thesis. RAY’I’RA corriputcs grovip limes, phase lime, 
signal losses, the cflccts of sporadic J; and elevation angles at both points. These data 
are saved for each successful propagation path and are passed to COiVlEFF. 

fhe raytracing algorithm computes the propagation path based data produced 
by NAIGEN. A ray trace w'hich ends within 1000 krn of the receiver site is saved for 
further processing. When two rays bracket the receiver, the program interpolates a ray 
that falls w'ithin some preset value, in this ca.se ten kilometers. The parameters describ¬ 
ing this ray arc saved for COMEFF. In the case wdiere only a single ray is found, the 
program again interpolates until a ray close to the receiver is found, fhe AMBCOM 
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User's Guide For Engineers [Ref. 15: pp.35-56] provides an in-depth treatment of 
raytrace algorithms. 

lUAYTRA estimates the amount of absorption for each ray. The absorption 
calculations are divided into four parameters as follows; 

• Lo is the divergence, i.e. free space spreading loss. 

• ionospheric absorption loss such as dcviative, nondeviative, and auroral losses. 

• L£j is the loss resulting from sporadic E. 

• Lq is the loss due to ground reflections along the path. 

The total loss for a path is the summation of all of these factors. Only ionospheric ab¬ 
sorption as it applies to auroral absorption will be described here. A full explanation 
of the RAYTIU'^ path loss algorithms may be found in Reference 15, pages 56 thru 80. 

2. Auroral Absorption. 

RAYTRA models the auroral absorption using a modified Foppiano model 
(Ref 15; pp. 64-71). The adaptations made to this model by SRI provide the model with 
some dependence upon the index. Auroral absorption is related to particle precipi¬ 
tation which is. in turn, heavily influenced by the activity of the geomagnetic field and 
A', is a descriptor of the current state of the planetary magnetic field. The sunspot 
number (which is used by the Foppiano model) does not correlate well with changes in 
the magnetic field and it is not a good indicator of the amount of auroral absorption 
present. 

The Foppiano model produces a median value for auroral absorption based 
upon the combined effects of substorm activity for a given hour of the day in a given 
month, and upon the level of solar activity. Solar activity is modeled with a 12 month 
running average of sunspot numbers which do not relate directly to geomagnetic field 
activity. Large variations in auroral absorption can be lost in the averaging process. 
On the other hand ihe addition of a large dependence upon A, could cause the model to 
produce incorrect results for quiet day predictions. SRI chose to include the A, index 
by the use of a normalized Gaussian distribution so that the location of the maximum 
amount of absorption is dependent upon but the magnitude of the maximum is not 
dependent upon A^. [Ref 15; pp. 64-65[ 

RA')TRA begins auroral absorption model by first calculating the factor Q, 
which is the one-way vertical absorption at 30 .VIHz that exceeds 1 dB. A riometer is 
used to measure changes in ionospheric absorption, at 30 MHz, in reference to quiet day 
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absorption levels. For RAYTRA the Q^ factor is calculated from tw'o other variables 
as 


Q] — Q\d+ Q\s' 


(13) 


= 21 d^x dr (14) 

and 

(15) 


where 

d/s^ = the drizzle/splash precipitation latitudinal dependence. 
dlsj= the drizzle,'splash precipitation time dependence, 
disft = the drizzle/spiash precipitation sunspot dependence, 
dlsg= drizzle,'splash precipitation longitudinal dependence and 
djs„ = the maximum latitude of drizzle,'splash precipitation. 

The term splash particle precipitation is associated with discrete ionospheric events that 
are limited spatially. Drizzle refers nearly constant and widespread particle precipitation 
that results from particle collisions in the upper atmosphere [Ref 22: p.322]. The sub¬ 
scripts Id and Is refer to drizzle and splash precipitation, respectively. The only pa¬ 
rameters modified for dependence are d^, s^. and dj, which are the latitude of the 
drizzle component, the latitude of the splash component, and the time dependence. The 
remainder of the Foppiano model has not been modified. [Ref 15 : pp. 64-66] 

The latitude of the drizzle component of the model is calculated with the 
equation [Ref 15: p. 66] 


d<t> = 


-L4>-4>n,f 

2al 


(16) 


where is the latitude of maximum absorption and is the latitudinal extent of the 
absorption zone. The A.’, influence is introduced with the equations [Ref 15: p. 66] 


4 > 


m 


67.864 - 0.6256/:^ ^Kp ^ 5 
75.616-2.176 ;Kp> 5 


(17) 


and 
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j 3.06 + .2765a; ;A; ^ 5 
[-.36+.96a; ;A'p>5 . 


( 18 ) 


The A^ inOuence for the time dependence is introduced similar to that for </)„. 
The time dependence is given by 


dj-= exp 


- [/, - - 
15.7 


(19) 


where 


9.9 - 0.46Ap ;A; > 5 

15.6-1.6a; ;Kp<S . 


( 20 ) 


Both and d-^ are normalized Gaussian functions. [Ref. 15: p. 67] 

The calculation of the splash component (s^) is also performed with a normal¬ 
ized Gaussian equation as follows: 


= exp 


2oi 



( 21 ) 


where 


66.799 - 0.9246A^ + (0.318 + 0.0828) | /1 \Kp^5 

78,256- 3.216Ap + (.28SA;,-.708)|t| \Kp> 5 ' 

_ r r-3 ;0< T< 15 
;15< r<24 


( 22 ) 


(23) 


and T is the corrected geomagnetic time. This component value along with the others 
is used to produce Q,. [Ref. 15: p.68[ 

is used to calculate the median absorption value y4„ , w^hich must be adjusted 
from a monthly median value to one corresponding to A,. The final value for absorption 
is obtained from 


where 




\ I. 


( 24 ) 
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}\.^= 1.8(0.17+ 0.14A». 


(25) 


The factor of 1.8 was included to scale the equation to match high latitude empirical 
data. [Ref 15: pp.69-71] 

3. Sporadic E Calculations 

R.‘\YTR.'\ offers two choices for computing sporadic E (£,). The fust method 
computes the reflection of signals for frequencies below the blanketing frequency 
The blanketing frequency is computed based upon the location of the control point. 
The second method performs E, reflection calculations for all frequencies less than f^E,. 
The value for f„E, is based upon the upper, median, or lower ITS/„£, coefficients passed 
from NATGEN. The choice of value is specified by the user. RAYTRA calculates/£, 
dependent upon the percentage of E, specified by the AMBCOM user. If 90% £, is 
specified then RAYTRA uses the upper decile value ofy^£, from the ITS file. For 50% 
£, RAYTRA uses the median decile value. None of these E, models were designed to 
simulated sporadic E at high latitudes. (Ref. 15: pp.53-56] 

RAYTR,‘\ offers two models for sporadic E (£,) in reference to absorption loss, 
the Sinno model and the Phillips model. Based upon the work of Hunsucker [Ref 7: p. 
351], the author decided to examine the effect of differing amounts of £, on the gener¬ 
ation of the model data. .The Sinno model was chosen because it allows the user to 
choose which of the decile values to use for f„E, for all sporadic E calculations. The 
Phillips model uses only the median decile value. [Ref 15: pp.71-77] 

4. RAYTRA Input Variables 

RAYTRA, like NATGEN, provides the user with a number of options for 
controlling program execution. Options are specified with numbered control cards in- 
stream following the program execution statement. Figure 8 on page 23 is an example 
of the RAYTRA controlling statements used for this thesis. 

Card two controls the amount of £, for a given execution. Three different types 
of input streams were chosen per campaign for comparing model accuracy as the per¬ 
centage of E, is varied. Figure 8 on page 23 is an example of a modified sporadic E in¬ 
put stream. The number 70 in columns 49-50 indicates a value of 70% f^Es will be used 
for all raytracing and loss calculations for sporadic E during the time period specified 
by the preceding NATGEN program. 
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Figure 8. RAYTRA input stream used to model winter site O propagation data. 

Project NONCENTRIC Site O data was compared to data generated by 
AMBCOM using 0''/o sporadic 0 calculations, 50% sporadic E, and modified sporadic 
E. For modified E„ the percentage was varied during the day. During the six hour pe¬ 
riod around midnight 90% E, was used, 70% during the three hour periods before and 
after the 90% window. RAY'fRA was allowed to default to 50% for all other hours 
during the day. For the 50% sporadic E runs, RAYTRA was allowed to default to the 
median ITS f,E, value. The sporadic E feature was turned off for the 0% E, input 
streams. 

Frequencies ‘'or RAYTRA are specified by control card four. Frequencies from 
Project NONCENTRIC were used in RAYTRA to enable automated comparison of 
model and measured data. 1 he automation of the comparison and statistical analysis 
process reduced the amount of time required to examine the results of a single execution 
and allowed more time for experimenting with the effects of changing control parame¬ 
ters. 
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E. COMEFF 


1. Overview 

The communications e/Tect (COMEFF) program evaluates the quality of a 
transmitted signal. Previously discussed programs, e.g., \ATGEN and RAYTRA, dealt 
with the ionosphere and possible ionospheric effects along a given propagation path. 
The only site-dependent information required for those two programs were the site lo¬ 
cations. CO.MEFF factors include transmit power, antenna configuration, signal band¬ 
width, and other communication site-dependent data. Based upon these local 
parameters and the data produced by RAYTRA, the performance of a particular com¬ 
munications link using given a frequency can be characterized in terms of the signal-to- 
noise ratio (SNR), field strengths and doppler spread. 

2 . Program Features 

COMEFF uses the data generated by RAYTRA for all of the saved modes as 
the basis for evaluation of a particular communications link. The raytrace descriptive 
data, for each ray, includes: 

• take-off angle, 

• group time. 

• phase time, 

• path loss, 

• noise power (per 1 Hz), 

• arrival angle and 

• transmitter frequency [Ref 15: p. 112]. 

Multiple raytraces may be described for each frequency. COMEFF combines the effects 
of multiple modes to produce a single reception statistic. 

COMEFF allows the user to include the effects of the antenna configuration 
on the communications link. The antenna data is taken from ANTLIB.DAT, the 
AMBCO.M antenna file. Unfortunately the antennas required for this thesis were not 
in the standard SRI antenna file. The antenna patterns for this thesis were taken from 
References 14 and 13 and adapted for use in AMBCOM. This is a manual process since 
AMBCO.M docs not provide the user with a program for entering the antenna data. 
Appendix A contains the antenna gain tables used for the summer and winter Clyde 
Rivcr-Leiccster campaigns. 

The contmunications effect program provides a variety of output options. 
CO.MEFF will calculate the signal-to-noise ratio (SNR), group time, doppler shift. 
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phase, and delay spread for all frequencies specified in RAYTRA. Other options provide 
field strengths or analysis of signal quality, e.g. the bit error rate. All requested data 
may be printed in several formats which include a listing of SNR, doppler shift, phase 
shift, and group times for each mode CO.MEFF also produces a single report that sums 
the data for each frequency at a particular time. 'I he SNR specified in the CO.MFFF 
summary report is a weighted accumulation of SNRs for all modes received from 
R,\^'TR,\. 

Weighted SNRs from the COMFIT summary report were compared to meas¬ 
ured SNR values. This composite SNR is calculated by computing the power for each 
ray. combining these values and then subtracting the noise value. The power A, for each 
ray is computed with the equation 


A, = ■ 


P,G,G,c^ 


4;rl0 


[T/'OJ 


f lO'^ lotC/ioj 


(26) 


u iiere 

p = transmitter power (W) 

O'. = transmitter antenna gain, 

O', = receiver antenna gain. 

L, = path loss for the ith ray, 

L. = receiving antenna loss and 
1' = Ircqucncy (.MI Iz). 

1 ho composite SNR \alue is given by the equation 

S.VR=101og,o(Xi'4,)-tVe, (27) 

where .V, is the noise power density in dBW. [Ref. 15 : p.83) 

3. COMEFF Input Variables 

Figure 9 on page 26 is an example of the COMEFF input stream. Two control 
cards are used to control program execution. The first card specifies transmit power, 
path loss threshold, and the use of antennas. Transmit power is given in units of watts 
per hertz. A bandwidth of 50 Hz and a transmitter power output of 400 watts results 
in 8 W Hz. The path loss threshold is used to exclude modes with a path loss greater 
than the specified threshold. This parameter was set at 400 dB which effectively allowed 
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all modes to be included. The second COM LIT control statement is used to specify 


antennas. 


$SET DEF DKAO;lAMBCOM] 

$ASSIGN SY's$INPUT FOR055: 

$ASSIGN DKAO:IWILSONlWAVEFORM.DAT FOn020: 
$ASSIGN DKAO:IWILSON)COMOUT.DAT FOn059: 
$ASSIGN DKAO:lAMBCOMlCOMEFF.DAT FOn007: 
$ASSIGN DKAO:tAMBCOMlANTLIB.DAT FOR006: 

$R DKAO:(AMBCOH)COMEFF 
NONCENTRIC 18 JANUARY lo''? 

0 0. 1600 1 i 0 8.0 

WHIP SLOPEVEC 


0 400 

1 
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in. PRO.IECT NONCENTRIC 


A. sn ES 

1. Overview 

Project NONCUN'lRiC was a study ol conwnunications at IIF frequencies at 
liigl) latitudes conducted by llic University of Leicester, U.K. I’hc concept of the project 
was to place a transniittcr in the polar cap area and to transmit IIF signals to receivers 
at several locations in polar cap and midlatitudes regions. Figure 10 shows approximate 
locations of transmitter and receiver sites, fable 1 on page 28 provides locations of sites 
in geographic latitude and longitude. 



1 he equipment configuration used for Project NONCEN'I RIC receiver sites 
varied from site to site No data arc available which accurately describes the receiving 
equipment, distribution systems, or antennas in operating environments. Data produced 
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Table 1. RECEIVER SITE LOCATIONS 


SITE 

LATITUDE 

LONGITUDE 

C 

57.10 N 

02.04 W 

D 

44.48 N 

68.46 W 

O 

53.00 N 

01.00 VV 


by AMBCOM were not adjusted to account for site-dependent anomalies. The data 
produced for sites C and D exhibited large average errors at higher frequencies, i.e. 30-40 
dB. It was assumed that site noise may have lowered the measured SNR, and the site 
C and D A.MBCOM SNRs were adjusted down 15 dB to account for noise in the re¬ 
ceiver system. 

The geometry of all the Project NONCENTRIC antennas are recorded, but 
none of the antenna patterns were measured at their specific site. No measurements 
were taken of site noise levels, ground permitivities, or ground conductivities. All an¬ 
tenna patterns and local noise levels used to generate AMBCOM data are estimated. 

2. Transmitter Site 

The transmitter site was located at Clyde River, Canada. Clyde River is located 
within the polar cap region at 70® north latitude and 69° west longitude. Propagation 
data used for this thesis are from midlatitude receiver stations so paths from Clyde River 
must pass through the auroral oval region. The Clyde River site transmitter power 
output level was 400 watts. 

Two different antenna configurations were used at the Clyde River site. During 
the summer campaign of 1988 a modified BUTTERNUT two-trap antenna (type 
HF6V-X) was used for all transmissions. The wdnter campaign employed a resistive 
whip antenna. These two antenna geometries were modeled by Tsolekas [Ref. 13; pp. 
11-16) using descriptive data provided by the University of Leicester. The antenna 
models were created using the Numerical Electromagnetics Code, Ver-ion 3 (NEC-3) 
using the Sommerfield-Norton approximation for a finite ground plane. The results from 
this modeling process were used to construct antenna gain tables used by AMBCOM. 

3. Receiver Sites 

, I 

Site O (University of Leicester) files were the only sources used for the earlier 
works by Tsolekas [Ref 13] and Gikas [Ref 14). In order to compare their results with 
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the data produced by A.MBCOM the same site data and antenna configurations were 
employed. This thesis also includes some comparisons of data generated by AMBCOM 
for sites C and D. 

Site O was equipped with a sloping vee antenna for all Project NONCENTRIC 
campaigns. This antenna was modeled by Gikas [Ref. 14: pp.10-16] using NEC-3. Gikas 
used the antenna geometry provided by the University of Leicester plus estimated values 
for permitivity and conductivity as inputs to NEC-3. The final form of the antenna gain 
tables for AMBCOM are interpolated from the results produced by Gikas for each of 
the Project NONCENTRIC frequencies. 

B. PROJECT NONCENTRIC DATA 

1. Signal Format 

Project NONCENTRIC data were taken from two of the four high latitude 
campaigns. The summer and winter campaigns were conducted 17 July-12 August 1988 
and 18 January -12 February 1989 respectively. The data were collected from three sites 
which are labeled as sites C, D, and O. 

Project NONCENTRIC signals were transmitted over a range of frequencies 
by stepping from 3.2 .MHz to 23.9 MHz at preset time intervals. These transmissions 
were formatted to provide a repetitive pattern of signal forms which could easily be re¬ 
cognized; Figure II on page 30 is a diagram of the signal format. The call sign (CS) 
was sent in morse code at five second intervals. The Barker code sequence consisted of 
a bipolar code with bit length of 250^5 repeated continuously during a 30 second inter¬ 
val. The Doppler signal represents 30 seconds of transmission with an unmodulated 
carrier. The call sign and Doppler signals were used for signal analysis by the automated 
signal recognition process used by the University of Leicester. 

2. Project NONCENTRIC Data Analysis 

An automated signal recognition program was used by the University of 
Leicester to validate the received transmissions and separated the data into the five files 
shown in Table 2 on page 30. The automated signal recognition process was validated 
by comparing the data from the five files with the results obtained by manual signal se¬ 
lection of a small subset of received transmissions. The automated signal recognition 
program was in agreement with manual selections about 80%-90% of the time [Ref 
23: p. 4], 
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Figure 11. Tlie time sequence of tlie Project NONCENI'RIC signals showing 
Doppler and call sign sequences. (From Ref. 23J 


Table 2. FROJECl MONCENI RIC WINTER/SITE O FILES 


Fll.li NAME 

FILE DFT IMTION 

OWINTER.SIG 

I’assed all test: valid data 

OWlNIER.rCS 

Passed only the call sign test 

OWINTER.r SI 

Passed only the spread index test 

OWINTER.INT 

Data judged to be interference 

OWIMER.NOS 

Data that failed all tests 


3. Spread Index Test 

The spread index test provides a simple method for judging the amount of 
Doppler spreading for each signal. 1 he method entails an analysis of the Doppler por¬ 
tion of the transmitted signal by observing the frequency shilt of an unmodulated carrier 
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with time. The amount of Doppler spreading is indicated the observed frequency shift 
of the unmodulated carrier during a 30 second transmission. A noise spectrum, assum¬ 
ing a uniform noise distribution, would not be expected to exhibit the same spectrum 
as an unmodulated carrier with Doppler spreading. The spread index test takes advan¬ 
tage of this property to separate possible valid signal data from noise data. [Ref 23: p. 
11 

The spread index test is a multiple step process. The received signal's spectrum 
is first obtained from a 1000 point Fast Fourier Transform (FFT). Figure 12 on page 
32 is an example of the resulting spectrum. The mean signal amplitude between -25 Hz 
and -12.5 Hz is taken as the mean noise level for the signal. This mean noise level value 
is subtracted from the spectrum to yield the form shown in Figure 13 on page 33. The 
area under the remaining signal is defined as the spread index. Since the spread index 
is proportional to the width of the remaining signal it is also proportional to the Doppler 
spreading of the signal. 

The spread index is used in conjunction with the signal strength to determine 
signal validity. The selection criteria are as follows: 

• A signal passes the spread index test if the spread index is large and the signal 
strength is greater than a preset threshold value. 

• The signal is invalid if the spread index is negative. 

• The signal is invalid if the spread index and signal strength are small. 

• The signal is valid if the spread index is small and the signal strength is large. 

All signals were passed to the call sign test regardless of the results of the spread index 
test. 

4. The Call Sign Test 

The call sign test provided a further means of validating data. The signal spec¬ 
trum was produced with a 1000 point FFT. The mean noise level (A',) was computed 
by examining the frequency spectrum from -50 Hz to -12.5 Hz. The value of the mean 
noise (A’j) used for the call sign test is increased to insure that a call sign signal mainlobe 
and sidelobes are present. The call sign mean noise value is given by 

A2 = 2.5AV (28) 

The value Aj is subtracted from the signal spectrum to produce the call sign signal. 
Figure 14 on page 34 is an example of the call sign speotrum after the mean noise level 
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Figure 12, The I’roject CONCENTRIC signal spectrum after the 1000 point FFT. 
(From Ref. 23] 

is removed. This signal is then examined lo determine if it matches the call sign pattern. 
(Ref. 23: p. 2] 

5, Intcrfeience 

The automatic signal recognition process was designed to detect interfering 
signals and to mark the transmission data accordingly. A signal is considered interfer¬ 
ence if it passes the spread index test but fails the call sign test. A signal-to-noise ratio 
(SNR) greater than 10 dB should pass the call sign test. If it fails the call sign test then 
the signal is considered to be interl'crence. A signal is considered invalid with a SNR less 
than 10 dB. [Ref 23: p.3| 

6. Signal Calegoriratiou 

1 he tested data arc categorized and stored in one of five files. 1 he file categor¬ 
ization criteria is as follows: 

• 1 he data passes both tests and is written to the signal (SIG) file. 
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Figure 13. The Project NONCLNTRIC FFT Spectrum «itli the average noise re¬ 
moved, IFrom Ref. 23] 

• I hc data failed both tests and is written to the nonsignal (NOS) file. 

• The data passes tlic spread index test onlv; the data is saved in the spread index file 
(ISI). 

I, 

• llie data passes the spread index test but is determined to be interference; the data 
is written to the interl'ercnce (INT) file. 

• The data passes tlic call sign test only; the data is written to the call sign (PCS) file. 

The data from the SIG files were used to compare with the AMBCOM generated data. 
(Ref 23: p. 3] 

Each of the receiver sites used for this thesis produced its own set of five data 
files, e.g.. SIG, NOS, TCS, FSl, and INI', fhe data are arranged in ascending order 
with respect to time. Table 3 on page 35 lists the attributes of each of the data elements. 

The PEAK and NOISE units are one of the major problems with the Project 
NONCEN I'RIC data. The two values (or these two elements are taken from.the 1000 
point EFT and they cannot be translated into the field strength and noise power values 















Figure 14. The Project NOINCENTRIC call sign spectrum produced by the 1000 
point FFT. (From Ref. 23) 

present at the antenna. Fortunately PEAK and NOISE are both in units of dB/^rV. 1 he 
signal SNR is obtained by subtracting NOISE from PEAK. T his SNR value was then 
compared to the AMBCOM SNR values produced for the COMEFF summary report. 

7. The Signals File Data 

T he end product of this collection and analysis was the SIG file, which repres¬ 
ents a very small portion of the total number of signals recorded during the Project 
NONCENTRIC campaigns. For example, the file OWIN'fER.SlG, the valid winter 
campaign signals received at site O, contains only 38 records out of a population of 609 
received signals at 3.2 MIIz. T he highest number of records for a particular frequency 
is 266 out of a population of 609 at a frequency of 14.4 MHz. The signal files used for 
this thesis are a small sample in comparison to the total number of signals that were 
received. 
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Table 3. PROJECT NONCENTRIC DATA ELEMENTS 


NAME 

Definition 

UNITS 

FREQ 

Frequency 

MHz 

HOUR 

Hour received 

UT 

MIN 

Minute received 

UT 

DAY 

Day of the week. 

none 

DATE 

Day of the month 


MONTH 

Month of the year 


YEAR 

Year 


PEAK 

Signal amplitude 

dB mV 

NOISE 

Mean noise level 

dB mV 

SPREAD 

Spread index 

dimensionless 


(I 
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IV. AMBCOM DATA ANALYSIS 


A. ANALYSIS METHODOLOGY 

1. General 

AMBCOM produces data which exhibits difTerent trends as frequency is 
changed. Sunimar}’ statistics, produced through manual analysis methods, tend to 
smooth over detail information about each frequency. An automated means of com¬ 
paring A.MBCOM data to measured data was designed to alleviate this problem. Com¬ 
puter programs compared A.MBCOM data to Project NOXCENTRIC data and 
accumulated statistics by frequency and sporadic E model. These statistics were then 
displayed graphically to highlight trends, by frequency, in AMBCOM's performance. 

2 . Data Synthesis 

Project XONCEXTRIC data is in a format that is easy to manipulate with a 
computer hut AMBCO.M data is not. A.MBCOM SNR data used for comparison was 
taken from the summary report produced by COMEFF. The data must first be removed 
from the report and stored in a form which simplifies the comparison process. Originally 
this task was performed manually by reviewing each line of the COMEFF summary re¬ 
port and erasing those elements that were not needed. A simpler method was pro¬ 
gramming the V.AX Workstation to perform this function which is one of the purposes 
of the modify data program (MODDAT). 

MODDAT is called inside .AMBCOM input streams which generate data for 
each campaign. It is executed once for each day or part of a day for which AMBCOM 
produces data. MODD.AT concatenates all COMEFF summary reports for that day 
and then removes the data required for comparison with the Project NOXCENTRIC 
data , e.g., the hour, frequency, and SNR. It also writes the month and day on each 
record. Files produced by MODDAT are allowed to accumulate until the end of 
A.MBCOM input stream execution, when they are all concatenated and sorted for later 
processing. 

3. Statistical Analysis 

The statistical analysis is a two step process. The first step compares the 
AMBCOM data with the actual (Project NOXCENTRIC) data. The matched records, 
actual records with no corresponding AMBCOM records, and AMBCOM records with 
no corresponding actual record are written to separate files for later processing. Statis- 
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tics regarding the average error, average SNR, number of matches, and number of misses 
of actual data are recorded. The final step accumulates statistics about the number of 
AMBCO.M records which did not match an actual record. 

The first step in the statistical analysis process is performed by the FORTRAN 
program \^'1LSTAT. WILSTAT matches A.MBCO.M data with actual data by com¬ 
paring the hour, month, day. and frequency of each record. Those AMBCOM records 
which do not match actual data are written to a file for later processing. Actual records 
which do not match AMBCO.M records are accumulated by frequency as misses. 

Each record that is matched is used to calculate, for each frequency, the average 
AMBCOM SNR and the number of matched records. All SNR values from A.MBCO.M 
matched data are accumulated in a table by frequency along with the total 

number of matched records for that frequency. At the end of the program this 

table was used to calculate the average SNR per frequency (SA7?„„) using the equation. 




X5.V 


XR 


"f'eq 


TOT, 


match 


(29) 


A file, ordered by frequency, is produced which lists the average SNR, total number of 
matched records, and the average error. 

The average error values are compiled from the error values recorded for each 
record. Once an A.VIBCOM record is matched to actual data, the error {ERR,} is cal¬ 
culated by subtracting actual SNR (S.\R„„) from A.VIBCOM SNR {SA7^„„(,„„). The 
frequency, month, day, hour, SA'R,„i,„„, error, and spread index of all matched 

records are written to a file for processing with a graphics routine. The average error 
per frequency (ERR,,,) is given by 


(S“«/) 

where TOT, \s the total number of records for that frequency. 

WILST.AT also calculates the percentage of matched records (Pror)- This value 
is the percentage of the total number of matched records (TOT) referenced to the total 
number of actual records (TOT,^,). The percentage is given by 
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p — TOT jQQ 

‘ TOT— TnT ^ 

‘ act _ 


(31) 


This equation was used by Gikas [Ref. 14 ] and Tsolekas [Ref 13] to produce their 
percentage of matched records. 

The second program of the statistical analysis process is called WILMAT., The 
purpose of this routine is to accumulate the number of unmatched AMBCOM generated 
predictions which were of acceptable transmission quality, e.g. AMBCOM predicts a 
strong signal for a time and frequency for which no signal was received. The key to this 
procedure is to select a proper threshold SNR value for acceptable signal quality. 
WILMAT performs this function by using the SNR„, values calculated by WILSTAT 
as the transmission quality threshold for each frequency. 

WILMAT calculates the percentage, by frequency, of matching AMBCOM data 
(Pf) referenced to the total number of AMBCOM predictions (TOT^) for a given fre¬ 
quency. All of the AMBCOM records used for this statistic must have a SNR greater 
than or equal to the SNR threshold for exceptable transmission quality. Pf is given by 



This percentage was included in the model analysis to provide some basis for judging the 
reliability of a given AMBCOM prediction. 

Using the data produced by WILSTAT, WILMAT produces the final summary 
report which lists by frequency 

• the average SNR, 

• the average error, 

• the total number of matched records, 

• the total number of unmatched actual records, 

• the total number of unmatched AMBCO.M records of acceptable transmission 
quality and 

• Pr 

Appendix C is the compilation of these statistics by site and campaign. All FORTRAN 
programs written to process the AMBCO.M data are listed in Appendix B. 
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B. AMBCOM INPUT STREAMS 

The 90°/o E, AMBCOVI input streams are designed to roughly emulate a pattern of 
£, frequency distributions that were described in Reference 7. The largest percentage 
of should occur at approximately midnight local time in the auroral oval, with less 
£, at later times. These conditions are very dependent on the location of the reflection 
points with respect to the auroral oval. is more intense within the auroral oval than 
at the lower latitudes. No attempt was made to manually adjust the percentage of E, 
at each of the reflection points. The percentage of £, was defined for the whole propa¬ 
gation path and A.MBCO.Vl adjusted for latitude and frequency differences. 

The purpose of the 501 o and 0°/o £, versions of the A.VIBCOM execution streams 
was to relate performance of the model with varying amounts of £, included in 
A.MBCOM calculations. The 50“'o E, data were modeled with using the median decile 
value for /,£,. The O^o sporadic E data were created without using A.MBCOM's £, 
model. 

AVIBCO.Vl £, models are based on midlatitude £, critical frequencies from the ITS 
file. The resulting £, models are not specifically designed to describe high latitude spo¬ 
radic E. The reason for this is that the original data produce by the Chatanika radar site, 
on which A.MBCO.M is based, did not have sufficient resolution to detect sporadic E 
(Ref 24j. The addition of a high-latitude £, model to A.MBCO.M is feasible and may 
be undertaken in the future; however, the effect of using the current AMBCOM sporadic 
E model for high latitudes had not been examined before this thesis. 

C. WINTER DATA 
1. Site O 

A.MBCO.M Site O data showed a decrease in percentage of matched records as 
the percentage of £, was decreased to 0. The 90° o, 50% and 0% E, data exhibited a 
PjoT of 88.9%, 80.5% and 68.3%, respectively. These percentages seem to indicate that 
the 90"o £, model is much more accurate in this panicular case, but this is not a com¬ 
plete statistical characterization. Figure 15 on page 40 shows the distribution of Pf for 
all three sporadic E models. Between 14.4 and 18.2 .MHz, the lower sporadic E models 
were able to match the actual data without predicting a large number of modes that had 
no Project NONCENTRIC equivalent. The 90% and 50% £, models showed a higher 
Pf at lower frequencies, i.e., 6.9 to 9.9 .MHz, than the data produced without sporadic 
E. The shapes of all three curves are approximately the same. 
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ire 15. The perceiit.ige of matching winter site O data relative to the total 
number of AM DCOM predictions (h) fi equency). 
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The average error also varied with the percentage of sporadic E. The largest 
amount of average error was exhibited by the 50% E, model, i.e. -9.6 dB. The 90% and 
O'/'o E, data average error was -7.6 and -7.0 dB, respectively. Figure 16 on page 42 shows 
that a large precentage of the error is distributed between -20 and + 20 dB, but the 50% 
sporadic E data shows a higher percentage of error outside of this range. Figure 17 on 
page 43 shows average error for all three models as it is distributed by frequency. The 
least amount of error is found in the frequency range between 9.9 and 20.9 .MHz. This 
is not a clear indicator of the performance of.AMBCOM since the accuracy of the an¬ 
tenna patterns used to produce the data is unknown. 

2. Site C 

Data generated by A.MBCO.M produced error patterns similar to site O. The 
percentage of matched records (Ptot) 91.7% (90% £,), 87.3% (50% EJ and 70.3% 
(0°'o £,). The distribution of P/is very different from the pattern of the site O data (see 
Figure 18 on page 44). P^ at higher frequencies did not greatly exceed those at lower 
frequencies and none of the percentages were larger than 55%. Some of these differ¬ 
ences may be caused by different receiving antennas at the two sites. 

The distribution of error (Figure 19 on page 45) exhibits that same skewed 
shape as site O data but the 50% E, curve does not have a second peak. The average 
error was -8.7 dB for 90% £„ -12.6 dB for 50% £, and -10.6 dB for 0% £,. The standard 
deviation of error for 90%, 50% and 0% £, models is 26.7 dB, 28.0 dB and 25.9 dB re¬ 
spectively. Figure 20 on page 46 shows the average error distributed by frequency. The 
average error is between -20 and +20 dB for frequencies between 6.8 and 23.2 MHz, 
roughly similar to the site O distribution. This is interesting since AMBCO.M is 
producing a smaller average error at the higher frequency inspite of a change in the re¬ 
ceive antenna. 

3. Site D 

The same trends continued for Site D data as were observed for Site O but the 
changes were not as dramatic. The overall percentage of matched records (£ror) does 
not vary' by more than a decade for all three modeling methods. The percentages are 
91.8%, 85.03o, and 81.8% for the 90%, 50% and 0% E, models respectively. The per¬ 
centage of matching A.MBCO.M records (see Figure 21 on page 48) showed a similar 
two peak pattern seen in site C data but peak values were in the 80%-90% range. The 
90 /0 £, model (Figure 21 on page 48) demonstrated the lowest percentage of matches 
per frequency of the three models. A close examination of the tabular statistics in Ap¬ 
pendix C indicates that the percentage of matches for 13.9 to 20.3 MHz improves with 
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Figure 16. The >vinfer site O aAerage error grouped in ten dB bins and plotted by 
percentage of error frequency. 
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Figure 17. Tlie nverage error, for winter site O, listed by frequeney. 
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18. The percentage of matching niiiler site C data relative to the total 
number of AMBCOM predictions (hy frequency). 
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Figiiie 19. The winter site C average eiror grouped in ten dB bins and plotted by 
percentage of error frequency. 
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Figuj e 20. I he average error, for winter site C, listed by frequency. 
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the 50® 0 and 0% sporadic E models relative to the 90% E, version. It is particularly 
interesting to compare site C and D data, because they both used the same type of an¬ 
tenna but produce very different results. This may be attributed in part to the location 
of two sites and the possibility that the two antennas do not perform the same at their 
respective locations. 

The average error for all three models displayed some variation. The 90% £, 
model produced an average error of -7.2 dB. The 50® o model exhibited an average error 
of -9.3 dB while the 0®/o model produced an average error of -11.1 dB. Figure 22 on 
page 49 shows the distribution of the average error for all three models. The curves are 
similar in shape and in location along the average error axis. Figure 23 on page 50 is 
the average error versus frequency for site D. Clearly, the site D distribution is similar 
to that of site C, and both sites demonstrate, as in the case of site O, that the smallest 
absolute value for average error is found at higher frequencies. 

D. SUMMER DATA 
1. Site O 

The summer campaign data produced results that varied slightly with the 
amount of sporadic E included in the AMBCOM model. The percentage of matched 
records (£ror) was, respectively, 95.0®/o, 84.8"/o and 79.1% for 90®/o, 50®/o and 0®/o £,. 
Figure 24 on page 51 is a plot of the percentage of matching A.MBCOM {Pf) data for 
each of the sporadic E models. AMBCOM predicted a larger percentage of the meas¬ 
ured data at high frequencies. The 90®/o, 50®/'o and 0®/o E, results are similar below 13.9 
MHz, but the 90®/o model exhibited lower Pf from 14.4 to 23.2 MHz. 

The average error of summer Site O data varied less than 3.5 dB between the 
three £, models. The 90®/o £, model produced an average error of -8.3 dB, but the 50®/o 
and 0®/'o models exhibited an average error of -11.05 and -11.87 dB, respectively. 
Figure 25 on page 53 shows the error distributions for all sporadic E models for the 
summer campaign. The distribution is skewed in the direction of negative average error. 
The standard deviations for these distributions are 18.2 dB (90®/o £, ), 19.3 dB (50®/o £,) 
and 23.0 dB (0®/o £,). The distribution of average error by frequency (Figure 26 on page 
54) is similar to the winter site O data. The lower frequency data exhibited larger ab¬ 
solute values of average error, but lower absolute error was the trend at higher fre¬ 
quencies, i.e., 13.9 to 20.9 MHz. There was very little variation in average error with 
£,. 
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Figure 21. The percentage of «inter site D niatcihng data r 
number of AMBCOM predictions (by frequency). 
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Figure 24. The perceiilage of matching summer site O data reference to the total 
number of AMBCOM predictions (In frequency). 
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Figure 25. Summer site O average error grouped in ten dB bins and plotted by 
percentage of error frequency. 
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ure 26. The average error, for summer site O, listed hy frequency. 


53 


AVERAGE ERROR (dB) 









2. Site C 


The summer site C data produced a high percentage of matched data {Pjot) for 
all three £, models. The respective values for 90%, 50% and 0% £, are 95.8%, 90.8” o 
and 81.8%. The distribution of was much different from that of the winter site C. 
Figure 27 on page 55 shows that Pf was smaller at higher I'requencies. This trend is the 
same for all three of the sporadic E models. For example, at 10.2 MHz all three of the 
models produced a high percentage of matching A.MBCOM records. Peaks and valleys 
in this figure roughly agree with a similar distribution for summer site D. The common 
antenna model may be a major contributing factor. 

Site C average error was obtained by subtracting 15 dB from the AMBCOM 
SNR to account for onsite noise at the receiving station. The average error is -6.2 dB, 
-8.9 dB and -13.3 dB for 90%, 50% and 0% £„ respectively. Figure 28 on page 56 is 
the error distribution for site C. The distribution is skewed along the negative error axis. 
The 0% distribution is offset from the 90% and 50% curves. The standard deviations 
for these error data are 15.9 dB, 16.4 dB and 21.5 dB for 90%, 50% and 0% £,. 

3. Site D 

Summer site D statistics were similar to those for summer site C. The percent¬ 
age of matched records were, respectively, 89.1%, 80.4% and 76.0% for the 90%, 50% 
and 0% £, models. The distribution of P; is plotted in Figure 30 on page 58. This dis¬ 
tribution exhibit a peak and valley pattern similar to that of summer site C. Values for 
Pf for site D are, generally, much higher than for site C. It is reasonable that the peaks 
and valleys are in part attributable to the behavior of the site antenna at different fre¬ 
quencies, i.e. the antenna patterns used for the AMBCOM modeling of these sites may 
not be completely correct. 

The average error for site D exhibited nearly the same pattern as for site C. The 
average error is -8.0 dB, -11.4 dB, and -17.0 dB for 90%, 50% and 0% £„ respectively. 
Figure 31 on page 59 shows the distribution of average error is centered near -20 dB for 
the 90% and 50% models. The 0% £, error distribution curve was shifted in the nega¬ 
tive direction along the average error axis with a peak value near -30 dB. The standard 
deviation of the error is 13.7 dB, 12.4 dB and 17.6 dB for the 90%, 50% and 0% models. 
Figure 32 on page 60 shows the distribution of average error by frequency. 
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Figure 27. The percent.nge of summer site C matching data reference to the total 
number of AMCCOM predictions fb} frequency). 
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Figure 28. Ihe suiiinier site C average error grouped in ten dB bins and plotted by 
percentage of error frequency. 
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Figure 29. Tlie average error, for suiiiiner si(e C, listed by frequency. 
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Figure 30. Tlic percentage of summer site D matching data reference to the total 
numher of AM DCOM predictions (by frcqueiicy). 
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Figure 31. l lie summer site 1 ) average error grouped in ten dB bins and plotted 


percentage of error frequency. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 

1. The Project NONCENTRIC Data 

Project NONCENTRIC data were not designed for the validation of 
computer-based propagation prediction models. For example, data did not contain 
necessarx' information such as 

• measured or computed antenna patterns, 

• measured received field strengths at the antenna, 

• measured signal and noise levels at the receiver, 

• ground constands at receiver sites, 

• measured receiver site noise levels and 

• periodic ionospheric data. 

For example, the antenna gain tables used by AMBCOM were modeled using estimated 
values of ground permitivity and conductivity for each receiver site. It is likely that these 
estimated parameters resulted in errors in the computation of AMBCOM SNRs. 
AMBCOM's validity as a high-latitude propagation prediction model cannot be fully 
established by Project NONCENTRIC data since the impact of these errors on 
AMBCO.M's accuracy cannot be evaluated. 

Project NONCENTRIC files were useful as a standard data base for comparing 
computer-based models. lONCAP, PROPHET and AMBCOM were executed using the 
same antenna gain tables, sunspot numbers, indices and frequencies. The data 
produced by these models were then compared to the same Project NONCENTRIC data 
using standardized analysis procedures. The results of these studies w’ere comparable 
since errors introduced by estimated input parameters were the same for each model. 

2. AMBCOM 

I he results from this thesis fully evaluate the usefullncss of AMBCOM for 
high-latitude propagation modeling, but AMBCOM does perform well when compared 
to lONC.'XP and PROPHET. As may be observed in Appendix D, the normalized me¬ 
dian error for AMBCOM, for the winter campaign, ranged from -.2 dB to 5.8 dB. The 
winter campaign normalized error for lON'CAP (Ref. 13; p. 24] and PROPHET [Ref 
14: p. 39] were 4.6 dB and .77 dB. respectively. AMBCO.M produced winter campaign 
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error distributions having standard deviations, depending upon £, model, ranging from 
15.9 dB to almost 20 dB. lONCAP's standard deviation for the winter campaign error 
distribution was 18.5 dB [Ref 13; p. 24j. PROPHET produced a standard deviation of 
14.7 dB for a similar distribution [Ref 14: p. 24|. Similar results were obtained for the 
site O summer campaign. 

A.MBCOM's batch input mode and operational flexibility make it ideal for au¬ 
tomated data analysis and communication link modeling with unique antenna config¬ 
urations. AMBCOM allow’s the user to 

• create large input streams representing several campaign days, 

• save data and report files for further processing, 

• add unique antenna patterns to A.MBCOM's antenna librar>', 

• include user written programs within the AMBCOM batch input stream and 

• review data produced bv program modules (e.g., review the NATGEN data passed 
to RAYTRA). 

lONCAP and PROPHET do not have all of these capabilities. lONCAP does not allow 
the user to easily add new antenna patterns to its antenna library'. Tsolekas [Ref 13] 
accounted for antenna effects by manually adjusting the lONCAP output. PROPHET 
does not provide a batch input mode, but PROPHET is an example of a interactive 
computer-based model with online editor features and both graphical and numerical 
output. Gikas [Ref 14] found it necessary to produce the PROPHET model data one 
hour and one frequency at a time for a 25 day campaign. 

A.MBCOM also has its limitations. AMBCO.M is difficult to use because it 
provides the user with little online support. Some of AMBCO.M's shortcomings noted 
during this thesis are 

• input control parameters must be located in specified columns, 

• there are no data entry' programs to assist the user in entering control parameters 
or antenna data, 

• there are no graphic output displays and 

• the programs run only on VAX VMS-based systems or a mainframe computer 
systems. 

Data produced by AMBCOM were displayed using a COMPAQ-386 and 
GRAPHTOOL because statistical trends are not easily visualized using numerical re¬ 
ports. 
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3. Statistics 

The accumulation of propagation prediction model statistics, by frequency, are 
useful for highlighting patterns in model performance. The analysis of lONCAP and 
PROPHET did not include statistics by frequency (Ref. 13 and 14]; however, the results 
from this thesis point to the importance of examining the performance of the model at 
each frequency. The general statistics and distributions of error provide only part of the 
picture of how well a model performs. For example, the 90% E, model for winter site 
O showed a Pjot nearly 90*lo and a low absolute value of average error. An exam¬ 
ination of the average error by frequency reveals that the absolute value of low fre¬ 
quency average error was very large. The distribution of Pf showed that a small 
percentage of A.MBCO.M predictions matched the measured data at low frequency. 
None of these trends may be observed with general statistics. 

B. RECOMMENDATIONS 

1. Comparison Data 

A new set of high-latitude data should be produced. This data should include 

• measured field strengths at the antenna and receiver, 

• measured noise levels throughout the receive site, 

• measured or, at least, carefully modeled antenna patterns and 

• periodic ionospheric data from the WWV broadcast and other solar-terrestrial 
sources. 

This will provide the model user with a standard set of transmission data and the back¬ 
ground information needed to create input for his computer-based model. 

2. AMBCOM 

A.MBCO.VI should be updated. The model should be redesigned with such fea¬ 
tures as 

• a menu-driven input editor, 

• a choice of operation modes, i.e., interactive or batch, 

• the capability to maintain a library of transmitter; receiver sites and 

• graphic output displays. 

One proposal is to use PROPHET'S user-friendly editor and graphic displays as a shell 
for the A.VIBCO.M software. The user would have the interactive features of PROPHET 
but would have access to the powerful A.MBCOM propagation prediction models. 
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AMBCOM should be downsized to run on PC-type microcomputers. The ma¬ 
jority of the militar}' user community has access to Zenith 248 type microcomputers. 
Placing AMBCO.M on a microcomputer will make it available to a larger segment of this 
user population. In addition, the microcomputer's graphics capability would provide 
AMBCO.M with a means of changing from numerical reports to graphical represent¬ 
ations of the data. The trends noted in this thesis were not dramatically noticeable until 
the data was graphed with PC-based software. AMBCOM provides no graphical dis¬ 
plays. 

3. Statistics 

All computer-based model examinations should use the same analysis method¬ 
ology. Statistics should include detailed information by frequency to include the average 
error and the percentage of matched records P, This will provide some insight regarding 
the behavior of the model at various frequencies. 
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Table 4 . RITSISTIVE WHIP GAIN TABLE 
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APPENDIX B. STATISTIC/.L ANAYLSIS PROGRAMS 


A. MODDAT 

^TfrTV*****^Tr**'A'*TV*******TVilt'!fr*-*Vr******Tfr*:fr***TiV:frVlr**Tfr**:fr**ifr*TVTVTV**yr'/f*TV****iVifrTVT*r*TV 

c* PROGRAM: MODDAT DATE: DEC 1990 

c* PROGRAMMER: D.J.WILSON 

c* MODDAT READS THE SUMMARY REPORT PRODUCED BY COMEFF AND STRIPS OUT THE 
c* THE TIME, FREQUENCY, AND SNR FOR LATER PROCESSING. THE DAY 
c* AND THE MONTH FOR WHICH THIS DATA WAS GENERATED IS 
c* ALSO RECORDED ON EACH OUTPUT RECORD. 

Q*TVTV**Vr*****Vf*‘A'^V*TV******iV******']V^V***yc**yrTV:fr**')ViV**yr*'^r******:fr:t****':VTViV**TV**:fryr 


REAL YFREQ,SNRIN ; SET UP REAL VARIABLES 

INTEGER MONTH,DAY,YHR ;SET UP INTEGER VARIABLES 

READ(58,2) MONTH,DAY ; READ IN THE MONTH AND DAY 

F0RMAT(I2,1X,I2) 

READ(5,10,ERR=300,END=500) YHR.YFREQ,SNRIN ;READ COMEFF LINE 
FORMAT(6X,I2,12X,F4. 1,49X,E11. A) ; IF ERROR GO TO 300 

IF (YFREQ. LT. 3. OR. YFREQ. GT. 24) GO TO 3 ; FREQ BETWEEN 3/24 MHz 

; NO - THEN GO TO 3 

IF (SNRIN.EQ.0) GO TO 3 ; YES - THEN IS SNR EQUAL TO 0 

; YES - THEN GO TO 3 

WRITE(4,20) YHR,YFREQ,SNRIN,MONTH,DAY ; NO - WRITE RECORD FOR 

; LATER PROCESSING 

F0RMAT(1X,I2,2X,F4. 1,2X,F7. 2,2X,I2,2X,I2) 


WRITE(4,20) YHR,YFREQ,SNRIN,MONTH,DAY 


GO TO 3 
GO TO 3 


;GO GET ANOTHER RECORD 

IN THE EVENT A FORMAT ERROR IS DETECTED DURING A READ 
THE PROGRAM TRANSFERS CONTROL TO STATEMENT 300 WHICH 
IN CAUSES THE PROGRAM TO GO BACK AND READ ANOTHER 
RECORD. THE ERROR RECORD IS EFFECTIVELY IGNORED. BY 
DOING THIS THE PROGRAM ONLY USES THOSE INPUT RECORDS 
WHICH MEET THE INPUT CRITERIA FOR A LINE OF INFORMATION 
FROM THE BODY OF THE COMEFF SUMMARY REPORT AND IGNORES 
THE REPORT HEADER AND TRAILERS. 


B. VVILSTAT 

c* PROGRAM: WILSTAT DATE: DEC 1990 

c* PROGRAMMER: D. J. WILSON 

c* THIS PROGRAM COMPARES THE AMBCOM DATA WITH THE PROJECT NONCENTRIC 
c* DATA. IF THE DATA MATCHES THEN THE SNR FOR BOTH RECORDS, ERROR, 
c* SPREAD INDEX, FREQUENCY, TIME, MONTH AND HOUR ARE WRITTEN TO 
c* FILE FOR LATER PROCESSING. THE FREQUENCY AND SNR FOR UNMATCHED AMBCOM 
c* DATA ARE WRITTEN TO A FILE FOR FURTHER PROCESS BY THE PROGRAM WILMAT. 
c* THE UNMATCHED NONCENTRIC RECORDS ARE COUNTED AND THE TOTAL IS 
c* MAINTAINED IN A TABLE BY FREQUENCY. ADDITIONALLY THE AVERAGE SNR BY 
C* FREQUENCY AND HOUR IS ACCUMULATED. 

Q****Vc**Vr*****Vr************** *********************************************** 

DATA SETUP *****iViV******'>V**->ViV*'iV*Vf**************-iV**ilr**')V** 
REAL*8 SNRTAB(16,24,3),ERRTAB(16,5),RECTOT(16) ; SET UP DOUBLE 

REAL*8 ACTEMP,MODTEMP,TOAVERR,TOTREC,PRCNT,TOTMOD ; PRECISION 

; VARIABLES. 
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REAL*4 MFREQ.MSNR ; SET UP SINGLE PRECISION 

; VARIALBLES. 

REAL*4 FREQ,SNR,SPI 

INTEGER MHR,MMON,MDAY,DOM,MON,XHR,HOUR ; IDENTIFY INTEGER VARIALBES 
ACTEHP=0. ; INITIALIZE 

('*■********* RFAD IN THF COMPARISON DATA ******ilr*5V*Vf****')V*********iir******"* 
READ(2:1^eJd=™o 5 FREQ:iSS,Spi:DOM,MON,XHR ; NONCENTRIC DATA 
10 F0RMAT(1X,F4. 1,IX,F7.2,1X,F7.2,IX,12,2X,11,2X,12) 

20 READ(3,30,END=1500) MHR,MFREQ,MSNR,MMON,MDAY ; AMBCOM DATA 

30 FORMATC1X,I2,2X,F4. 1,2X,F7. 2,2X,I2,2X,12) 


40 


MSNR=MSNR - 15.0 

M0DTEMP=0. 

L=1 

1=1 

J=1 

IFCMFREQ. EQ. SNRTAB(I,J,L)) THEN 
GO TO 100 

ELSE 


USED TO ADD 15 dB OF NOISE TO 
SITES C AND D DATA. 

INITIALIZE TEMP SPACE 
INITIAILIZE INDICES 


IS THE AMBCOM FREQ IN THE TABLE? 
YES - GO PROCESS 
NO - FIND AN EMPTY SPOT 
IS THIS AN EMPTY LINE? 

YES - SAVE FREQ 
GO PROCESS 
NO - 

INCREMENT BY 1 

GO LOOK AT NEXT LINE OF TABLE 


IF(SNRTAB(I,J,L). EQ. 0) THEN 
SNRTAB(I,J,L)=MFREQ 
GO TO 100 

ELSE 

1 = 1+1 
GO TO 40 
END IF 

END IF 

C ***** add SNR TO ACCUM SNR BY FREQ ************************************ 
100 J=MHR+1 

SNRTAB(I,J,2)=MSNR + SNRTAB(I,J,2) 

SNRTAB(I,J,3)=SNRTAB(I,J,3)+1 

c****** MATCH THE RECORDS ********************************************** 
105 ACTEHP=(M0N*1E5)+(D0M* 1E2) + XHR 
MODTEMP= (MM0N*1E5)+(MDAY*1E2)+MHR 
no IF(ACTEMP-MODTEMP) 200,150,300 ;DO THE MONTHS, DAYS AND HOURS 
150 1=1 ; MATCH? 

IF(FREQ-MFREQ) 200,155,300 ;YES - DO THE FREQ MATCH? 

155 ERR=MSNR-SNR ;YES - CALCULATE ERROR AND PRINT 

WRITE(4,160) MFREQ,MSNR,SNR,ERR,SPI,MON,MDAY,MHR 
160 FORMATCIX,F4.1,1X,F7.2,1X,F7.2,1X,F7.2,1X,F7. 2,IX,12,IX,12,IX,12) 

C******* INCREMENT THE ERROR TABLE AT THE CORRECT FREQ ****************** 
170 IFCMFREQ. EQ. ERRTABC1,1)) THEN ;DOES THE AMBCOM FREQ = TABLE FREQ? 
GO TO 180 ; YES - GO PROCESS 

ELSE ; NO - 

IFCERRTABCI,1).EQ. 0) THEN ;IS THIS AN EMPTY SLOT IN TABLE? 

ERRTABCI,1)=MFREQ ;YES - SAVE FREQ AND GO PROCESS 


180 


ELSE 

1 = 1+1 
GO TO 170 

ENDIF 

ENDIF 

ERRTAB C1,2)=ERRTAB C1,2)+ERR 

ERRTABCI,3)=ERRTABCI,3)+1 
ERRTABCI,4)=ERRTABCI,4)+ MSNR 


,NO - INCREMENT AND GO TRY AGAIN 


ADD ERROR FOR AVERAGE ERROR 

CALCULATION LATER 

ADD ONE TO MATCH RECORD COUNT 

ADD SNR FOR AVERAGE SNR CALCULATION 
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n n o n n 


READ(2,10,END=1400) FREQ,SNR,SPI,DOM,MON,XHR ;GET MORE ACTUAL DATA 
GO TO 20 ;GO GET MORE AMBCOM DATA 

C **** NO MATCH - GET ANOTHER ACTUAL DATA RECORD ****************** 


200 

210 


1=1 

IF(FREQ.EQ. ERRTAB(I,1)) THEN 
GO TO 250 

ELSE 


IF(ERRTAB(I,1).EQ. 0) THEN 
ERRTAB(I,1)=FREQ 
GO TO 250 


LIST MISS FOR ACTUAL RECORD IN TABLE 
DOES ACTUAL FREQ EQ TABLE? 

YES - GO PROCESS 
;N0 - 

;IS THIS SLOT EMPTY? 

; YES - SAVE ACTUAL FREQ IN TABLE 
;G0 PROCESS 


ELSE 


1=1+1 ;INCREMENT BY ONE GO GET MORE 

GO TO 210 ;DATA 

END IF 


END IF 

250 ERRTAB(I,5)=ERRTAB(I,5)+1 ; ADD ONE TO ACTUAL MISSES 

WRITE(9,10) FREQ,SNR,SPI,DOM,MON,XHR ; FOR THAT FREQUENCY 
READ(2,10,END=I400) FREQ,SNR,SPI,DOM,MON,XHR ; GET ANOTHER ACTUAL 
ACTEMP=0 ;RECORD AND GO AGAIN 

GO TO 105 

C ***** NO MATCH - GET MORE MORE MODEL DATA **************************** 
300 WRITE(7,1450) MFREQ,MSNR ;SAVE AMBCOM DATA FOR WILMAT 
360 GO TO 20 ;GET ANOTHER AMBCOM RECORD 

C 

C*** OUT OF ACTUAL DATA READ ALL REMAINING MODEL RECORDS TO NO MATCH TABLE 
1400 WRITE(7,1450) MFREQ.MSNR 
1450 F0RMAT(1X,F4.1,2X,F7.2) 

1460 READ(3,30,END=2000) MHR,MFREQ,MSNR,MON,NDAY 
C MSNR=MSNR-15.0 

GO TO 1400 


C 


C*** OUT OF MODEL DATA READ ALL REMAINING ACTUAL RECORDS TO NO MATCH TABLE 
1500 1=1 

1510 IF (FREQ. EQ.ERRTABC1,1)) THEN ;IS THIS ACTUAL FREQ THE CURRENT 
GO TO 1550 ;TABLE FREQ - YES THEN GO PROCESS 

ELSE 


IF(ERRTAB(I,1).EQ.0) THEN ;NO - IS THIS AN EMPTY SPOT? ■ 
ERRTAB(I,1)=FREQ ;YES - SAVE THE ACTUAL FREQ " 

ELSE 


1=1+1 ;NO - INCREMENT BY ONE AND GO 

GO TO 1510 ;LOOK AGAIN. 

ENDIF 


ENDIF 

1550 ERRTAB(I,5)=ERRTAB(I,5)+1 ; INCREMENT ACTUAL MISS COUNT FOR 

WRITE(9,10) FREQ,SNR,SPI,DOM,MON,XHR ;THIS FREQ THEN WRITE THE 
READ(2,10,END=2000) FREQ,SNR,SPI DOM,MON,XHR ; ACTUAL RECORD AND 
GO TO 1500 ;READ A NEW ACTUAL RECORD 


**+ THE END ROUTINE - SPIT OUT ALL THE ACCUM DATA FROM THE TABLES 


2000 DO 2100 1=1,16 

IF(ERRTAB(I,1).EQ.0) THEN 
CONTINUE 
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ELSE 

IF (ERRTAB(I,3). GT. 0) THEN 

AVEERR=(ERRTABC1,2)/ERRTAB(1,3)) 

AVESNR1=(ERRTABC1,4)/ERRTABC1,3)) 

TOAVERR= CAVEERR * ERRTABCI,3))+TOAVERR 
TOTMOD = ERRTABCI,3)+T0TM0D 
TOTREC= TOTREC + ERRTABC1,5) 

WRITEC5,2010) ERRTABC1,1),ERRTABC1,3),AVEERR,AVESNR1, 
C ERRTABCI,5) 

2010 F0RMATC1X,F4. 1,2X,F7.0,2X,F7. 2,2X,F10. 2,2X,F7. 0) 

ELSE 

AVERR=0.0 
AVESNR1=0.0 

TOTREC = TOTREC + ERRTABC1,5) 

WRITEC5,2010) ERRTABCI,1),ERRTABCI,3),AVEERR,AVESNR1, 
C ERRTABCI,5) 

ENDIF 
END IF 

DO 2170 J =1,24 
IF CSNRTABCI,J,3).GT. 0) THEN 

SNRTABCI,J,2)= CSNRTABCI,J,2)/SNRTABCI,J,3)) 
RECTOTCI)= RECTOTCD + SNRTABC I, J, 3) 

ELSE 

SNRTABCI,J,2)= -999. 

ENDIF 

2170 CONTINUE 
2100 CONTINUE 
2200 DO 2500 1=1,16 

IF CSNRTABCI,1,1).GT.0) THEN 

WRITEC8,2040) SNRTABCI,1,1),CSNRTABCI,J,2), J=l,12) 

2040 FORMATC 1',F4.1,12C1X,F4.0)) 

WRITEC 8,2050) SNRTABC1,1,1),CSNRTABCI,J,3),J=1,12) 

2050 FORMATC'2',F4.1,12C1X,F4.0)) 

WRITEC8,2042) SNRTABCI,1,1).(SNRTABCI,J,2), J=13,24) 

2042 FORMATC'3' ,F4'. 1,12C 1X,F4. 0)) 

WRITEC8 2052) SNRTABCI,1,1),(SNRTABCI,J,3),J=13,24) 

2052 FORMATC 4',F4.1,12C1X,F4.0)) 

WRITEC 8,2045) SNRTABC1,1,1),RECTOTCI) 

2045 FORMATC 5',F4. 1,IX,'TOTAL NUMBER OF RECORDS',2X,F7.2) 
ELSE 

CONTINUE 

ENDIF 

2500 CONTINUE 

TOAVERR= CTOAVERR/TOTMOD) 

PRCNT= CTOTMOD/CTOTMOD+TOTREC))*100 
WRITEC6,2600) FRONT, TOAVERR 
2600 FORMATCIX,F5.2,IX,F7. 2) 

END 


C. WILMAT 

c* PROGRAM: WILMAT DATE: DEC 1990 

c* PROGRAMMER:D.J. WILSON 

c* THIS PROGRAM EXAMINES THE UNMATCHED AMBCOM PREDICTIONS AND DETERMINES 
c* IF THEY ARE SIGNAL QUALITY TRANSMISSIONS. IF THE AMBCOMB PREDICTION 
c* PASSES THIS TEST, THEN THE COUNTER FOR THE TOTAL NUMBER OF AMBCOM 
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c* PREDICTIONS FOR THAT FREQUENCY IS INCREASED BY ONE. ONCE ALL THE 
c* AMBCOM PREDICTIONS ARE PROCESSED A REPORT CONTAINING, BY FREQUENCY, 
c* AVERAGE SNR, AVERAGE ERROR, NUMBER OF MATCHED RECORDS, NUMBER OF 
c* UNMATCHED ACTUAL RECORDS, NUMBER OF UNMATCHED AMBCOM RECORDS AND THE 
C* PERCENTAGE OF MATCHED AMBCOM RECORDS REFERENCED TO THE TOTAL NUMBER 
c* OF AMBCOM RECORDS FOR A GIVEN FREQUENCY. 


REAL*8 SNRTAB(16,6) ; SET UP DOUBLE PRECISION TABLE 

REAL*4 MFREQ,NOSNR,NOMAT ; RESERVE SPACE FOR REAL VARIABLES 

1=1 ; INITIALIZE INDICES 

J=1 

C********** BUILD UP TABLE OF AVERAGE FREQUENCY DATA *********yryr** 
10 READ(3,50,END=100) FREQ.NOMAT,AVEERR,AVESNR1,ACTNOMAT ; READ 


50 F0RMAT(1X,F4. 1,2X,F7. 0,2X,F7. 2 
60 IF(FREQ.EQ.SNRTABC1,1)) THEN ; 

GO TO 10 ; 

ELSE 

IF(SNRTAB(I,1).EQ.0) THEn’ 
SNRTABCI,1)=FREQ ; 
SNRTABCI,2)=AVESNR1 ; 

SNRTABCI,3)=AVEERR ; 
SNRTABCI,4)=NOMAT ; 
SNRTABCI,5)=ACTNOMAT ; 
GO TO 10 ; 

ELSE ; 

1 = 1+1 ; 
GO TO 60 ; 

ENDIF 


AVERAGE ERROR, AVERAGE SNR, 

NUMBER OF MATCHES, NUMBER OF 
ACTUAL MISSES. 

,2X,F10. 2,2X,F7. 0) 

IS THIS FREQ IN THE TABLE? 

YES - GO PROCESS 
NO - FIND AN EMPTY SPACE 
; IS THIS SPACE EMPTY? 

YES - MOVE FREQ TO FIRST SPOT OF 
THIS LINE OF THE TABLE AND THEN 
MOVE THE AVERAGE SNR, AVERAGE ERROR, 
NUMBER OF MATCHES PER FREQ, AND THE 
NUMBER OF UNMATCHE NONCENTRIC RECS 
GO PROCESS 

.NO - THIS SPACE IS NOT EMPTY 
INCREMENT INDEX BY 1 
GO LOOK AGAIN 


ENDIF 

C*.v***,v USE TABLE OF AVERAGE DATA TO SIFT UNMATCHED AMBCOM DATA ******* 


100 

150 


200 


300 


350 


READC4,150,END=380) MFREQ.NOSNR 
FORMATC1X,F4. 1,2X,F7.2) 

1 = 1 
J=1 

IFCSNRTABCI,1).EQ,MFREQ) THEN 


GO TO 300 

ELSE 

IFCI.EQ. 16) THEN 
GO TO 100 

ELSE 

1 = 1 + 1 
GO TO 200 

ENDIF 

ENDIF 

IFCSNRTABCI,2).GT.NOSNR) THEN 
GO TO 100 

ELSE 

CONTINUE 

ENDIF 

S.NRTABC I ,o)= SNRTABC 1,6) + l 


READ AN UNMATCHED AMBCOM RECORD 

INITIALIZE THE INDICES 

DOES THE UNMATCHED AMBCOM FREQ 
EQUAL THE FREQ IN THIS LINE 
OF THE AVERAGE DATA TABLE? 

YES - GO PROCESS 
NO 

HAVE WE EXHAUSTED THE TABLE? 

YES - GO GET ANOTHER RECORD 
NC - INCREMENT THE INDEX AND 
SEARCH NEXT LINE OF THE TABLE. 


IS AVERAGE SNR GREATER THAN 
UNMATCHED AMBCOM RECORD SNR? 
YES - THROWAWAY, GET A NEW REC. 
NO - GO PROCESS 


ADD ONE TO THE COUNT OF UNMATCHED 
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; AMBCOM PREDICTIONS WHICH HAVE NO 
; CORRESPONDING NONCENTRIC REC. 

GO TO 100 ; GET ANOTHER UNMATCHED AMBCOM REC. 

C **** OUTPUT THE TOTAL MISSES *******iv*vr->v******vr*****yr-jv***iv****5v****iv** 
380 DO 400 1=1,16 ; WRITE THE TABLE WHICH CONTAINS FOR EACH 

; FREQUENCY THE AVERAGE SNR, AVERAGE ERROR, 

; NUMBER OF MATCHED RECORDS, NUMBER OF 
; NONCENTRIC RECORDS WITHOUT MATCHES, AND 
; NUMBER UNMATCHED AMBCOM RECORDS. 

IF (SNRTAB(I,1).GT.0) THEN 

WRITE(5,2030) SNRTAB(I,1),SNRTAB(I,2),SNRTAB(I,3), 

C SNRTAB(I,4),SNRTAB(I,5),SNRTAB(I,6) 

2030 F0RMAT(1X,F4. 1,2X,F7. 2,2X,F7. 2,2X,F7.0,2X,F7. 0, 

C 2X,F10.2) 

ELSE 

CONTINUE 

END IF 

400 CONTINUE 
END 
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APPENDIX C. TABLES OF SITE STATISTICS 


Table 7. WINTER SITE O - 90% ES 


FREQ 

SNR 

ERR 

NUMBER 

ACTEAE 

A.MBCOM 


(MHz) 

(dB) 

(dB) 

MATCHED 

.MISSES 

MISSES 


3.2 

-63.4 

-83.2 

38.0 

0.0 

187.00 

16.9 

4.5 

-27.4 

-57.3 

17.0 

4.0 

202.00 

07.6 

6,8 


-28.6 

208.0 

0.0 

134.00 

60.8 

6.9 

3.3 

-17.8 

167.0 

0.0 

120.00 

58.2 

9.9 

18.0 

mmiii 

149.0 

0.0 

183.00 

44.9 

10.2 

sm 

1.3 

167.0 

1.0 

209.00 

44.4 

13.9 

24.3 

0.0 

200.0 

28.0 

140.00 

58.8 

14.4 

26.5 

2.2 

238.0 

28.0 

105,00 

69.3 

17.5 

28.1 

-0.6 

169.0 

57.0 

47.00 

78.2 

18.2 

29.1 

0.4 

170.0 

56.0 

36.00 

82.5 

20.3 

26.8 

0.2 

119.0 

30.0 

66.00 

64.3 

20.9 

24.3 

-1.9 

96.0 

11.0 

86.00 

52.7 


16.4 

-16.6 

22.0 

4.0 

81.00 

21.4 


PERCEM AGE MATCHED 88.93% AVERAGE ERROR -7.6dB 





































































































Table 8. WINTER SITE O - 50% ES 


FREQ 

SNR 

ERR 

NUMBER 

ACTUAL 

AMBCOM 

°/o 

(MHz) 

(dB) 

(dB) 

MATCHED 

.MISSES 

MISSES 


3.2 

-61.4 

-81.2 

38.0 

0.0 

167.00 

18.5 

4.5 

-28.0 

-57.8 

17.0 

4.0 

198.00 

7.9 

6.8 

-7.5 

-32.5 

208.0 

0.0 

167.00 

55.5 

6.9 

0.7 

-20.4 

167.0 

0.0 

118.00 

58.6 

9.9 

11.7 

-10.0 

141.0 

8.0 

167.00 

m 

10.2 

11.9 

-8.5 

157.0 

11.0 

228.00 

40.8 

13.9 

28.8 

4.1 

144.0 

84.0 

79.00 

64.6 

14.4 

31.4 

6.8 

176.0 

90.0 

59.00 

74.9 

17.5 

30.2 

1.5 

150.0 

76.0 

25.00 

85.7 

ISH 

30.4 

1.7 

160.0 

66.0 

29.00 

84.7 

20.3 


1 

o 

118.0 

31.0 

65.00 

64.5 

20.9 

24.0 

-2.3 

96.0 

11.0 

81.00 

54.2 


16.4 

-16.6 

22.0 

4.0 

76.00 


PERCENTAGE .MATCHED 80.55°/o AVERAGE ERROR -9.61 dB 
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Table 9. WINTER SITE O - 0% ES 


FREQ 

SNR 

ERR 

NUMBER 

ACTUAL 

AM BCDM 

".0 

(MHz) 

(dB) 

(dB) 

MATCHED 

MISSES 

MISSES 


3.2 

-67.9 

-88.0 

37.0 

1.0 

178.00 

17.2 

4.5 

-44.0 

-74.5 

14.0 

7.0 

224.00 

5.9 

6.8 

-5.5 

-31.1 

132.0 

76.0 

100.00 

56.9 

6.9 

0.1 

-21.6 

lOS.O 

59.0 

109,00 

49.8 

9.9 

17.8 

-3.2 

96.0 

.-3.0 

161.00 

37.4 

10.2 


-1.1 

93.0 

75.0 

186.00 

33.3 

13.9 

29.1 

4.2 

144.0 

84.0 

83.00 

63.4 

14.4 

31.6 

6.9 

1 75.0 

91.0 

58.00 

75.1 

17.5 

30.0 

1.4 

152.0 

74.0 

29.00 

84.0 

18.2 

30.2 

1.5 

164.0 

62.0 

29.00 

85.0 

B^m 

26.4 

-0.3 

IIS.O 

31.0 

63.00 

65.2 

20.9 

23.8 

-2.4 

97.0 

10.0 

83.(.X) 

53.9 

23.2 

16.6 

-16.4 

22.0 

4.0 

77.00 

22.2 

PERCENTAGE MATCHED 68.32% AVERAGE ERROR -7.04 dB 
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Table 10. WINTER SITE C - 90% ES 


FREQ 

SNR 

ERR 

N EMBER 

ACTUAL 

AMBCOM 

°/o 

(VI Hz) 

(dB) 

(dB) 

MATCHED 

MISSES 

MISSES 


3.2 

-54.6 

-73.0 

22.0 

0.0 

234.00 

8.6 

4.5 

-14.8 

-36.6 

19.0 

3.0 

180.00 

9.5 

0.8 

8.7 

-12.5 

138.0 

0.0 

133.00 

50.9 

6.9 

13.5 

-5.3 

161.0 

0.0 

139.00 

53.7 

9.9 

17.7 

-2.5 

80.0 

1.0 

259.00 

23.6 

10.2 


2.7 

98.0 

2.0 

178.00 

35.5 

13.9 

26.9 

5.7 

43.0 

4.0 

194.00 

18.1 

14.4 


3.9 

57.0 

7.0 

290.00 

16.4 

17.5 

WSBM 

-11.3 

90.0 

23.0 

254.00 

26.2 


7.9 

-17.3 

118.0 

34.0 

223.00 

34.6 

20.3 

15.7 

-5.8 

15.0 

2.0 

199.00 

7.0 

20.9 

32.9 

13.0 

6.0 

0.0 


4.9 

23.2 

12.7 

-1.9 

1.0 

1.0 

169.00 

0.6 

PERCENTAGE MATCHED 91.683'o AVER.^GE ERROR -8.75 dB 
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Table 11. WINTER SITE C - 50% ES 


FREQ 

SNR 

ERR 

NUMBER 

ACTUAL 

A.MBCO.M 

o/o 

(MHz) 

(dB) 

(dB) 

MATCHED 

MISSES 

MISSES 


3.2 

-55.5 

-71.8 

11.0 

11.0 

85.00 

11.5 

4,5 

-20.4 

-41.1 

7.0 

15.0 

88.00 

7.4 

6.8 

5.6 

-14.8 

58.0 

80.0 

63.00 

47.9 

6.9 

8.8 

-9.2 

72.0 

89.0 

87,00 

45.3 

w^m 

3.7 

-16.5 

27.0 

54.0 

226.00 

10.7 

1U.2 

8.2 

-10.5 

30.0 

70.0 

189.00 

13.7 

13,9 

25,7 

n 

16.0 

31.0 

39,00 

29.1 

14.4 

23.7 

3.3 

13.0 

51.0 

61.00 

17.6 

17.5 

9.6 

-13.6 

35.0 

78.0 

59.00 

37.2 

18.2 

10.6 

■-14.7 

54.0 

98.0 

43.00 

55.7 

20.3 

11.3 

-11.3 

7.0 

10.0 

72.00 

8.9 

Eom 

15.7 

-3.6 

2.0 

4.0 

64.00 

3.0 

Eom 

0.00 

-3.57 

0.0 

2.0 

89.00 

0.0 

PERCENTAGE MATCHED 35.89% AVER.^GE ERROR -13.87 dB 
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Table 12. WINTER SITE C - 0% ES 


EREQ 

SNR 

ERR 

NUMBER 

ACTUAL 

AMBCOM 


(MHz) 

(dB) 

(dB) 

MATCHED 

MISSES 

MISSES 


3.2 

-61.7 

-S0.4 

21.0 

1.0 

191.00 

9.9 

4.5 

-25.0 

-47.6 

17.0 

5.0 

134.00 

11.3 

IBHI 

5.5 

-15.6 

86.0 

52.0 

73.00 

54.1 

6.9 

2.7 

-15.6 

105.0 

56.0 

161.00 

39.5 

9.9 

13.9 

-6.7 

56.0 

25.0 

218.00 

20.4 

10.2 

21.0 

1.0 

59.0 

41.0 

134.00 

30.6 

mmm 

27.6 

6.6 

32.0 

15.0 

86.00 

27.1 

14.4 

25.2 

5.7 

50.0 

14.0 

119,00 

ED 

tarn 

15.7 

-7.S 

86.0 

27.0 

103.00 

45.5 

IS,2 

17.7 

-7.8 

116.0 

36.0 

86.00 

57.4 

20.3 

15.8 

-5.7 

15.0 

2.0 

151.00 

9.0 

20.9 

mm 

13.0 

6.0 

0.0 

121,00 

4.7 

23.2 

12.7 

-1.9 

1.0 

1.0 

175.00 

0.6 


PERCENTAGE MATCHED 70.27% AVER.^GE ERROR -10.59 dB 
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Table 13. WINTER SITE D - 90% ES 


FREQ 


(Mil/) 


3.2 


.5 


6.S 


SNR 

ERR 

(dB) 

(dB) 


-57.1 

-17.8 

-36.6 

9.2 

-16.0 

14.1 

-8.3 

17.7 

-4.9 

22.7 

0.1 

IS.3 

-7.3 

19.6 

-7.4 

17.7 

-14.6 


3.0 


107.0 


71.0 


195.0 


206.0 


250.0 


220.0 


1.0 


250.0 


286.0 


192.0 


197.0 


172.0 


percentage: matched 91.81% average error -7.20 dB 



ACTEAL 

AMBCO.M 

8o 

MISSES 

MISSES 


0.0 

160.00 

1.8 

30.0 

255.00 

29.6 

0.0 

65.00 

85.1 

0.0 

184.00 

51.5 

5.0 

204.00 

50.2 

2.0 

128.00 

66.1 

27.0 

199.00 

52.5 

42.0 

194.00 

54.4 

44.0 

185.00 

57.5 

31.0 

173.00 

62.3 

27.0 

109.00 

63.8 

14.0 

72.00 

73.2 

17.0 

74.00 

69.9 
































































































Table 14. WINTER SITE D - 50% ES 


FREQ 

SNR 

ERR 

M MBER 

ACTUAL 

A M BCDM 

1 

‘(o 

(.MHz) 

(dB) 

(dB) 

MATCHED 

MISSES 

.MISSES 


3.2 

-37.9 

-58.8 

3.0 

0.0 

152.00 

1.9 

4.5 

-18.8 

-37.6 

107.0 

30.0 

242.00 

30.7 

6.8 

wsm 

-18.4 

368.0 

3.0 

61.00 

85.8 

6.9 

11.3 

-11.1 

194.0 

1.0 

162.00 

54.5 

9.9 

12.3 

msm 

174.0 

37.0 

177.00 

49.6 

10.2 

14.6 

-7.3 

190.0 

62.0 

144.00 

56.9 


18.3 

-7.0 

199.0 

48.0 

86.00 

69.8 

14.4 

19.7 

-7.2 

214.0 

59.0 

94.00 

69.5 

17,5 

16.2 

_1 

-16.6 

228.0 

66.0 

51.00 

81.7 

18.2 


-12,7 

254.0 

63.0 

50.00 

83.6 

20.3 

34.0 

3,4 

188.0 

31.0 

74.00 

71.8 

20.9 

37.5 

5.9 

192.0 

19.0 

69.00 

73.6 


43.0 

8.6 

169.0 

20.0 

73.00 

69.8 


PERC'EM AGH MATCHED 84.%'% AVERAGE ERROR -9.26 dB 




























































































Table 15. WINTER SITE D - 0% ES 


EREQ 

SNR 

ERR 

NUMBER 

ACTUAL 

A.MBCOM 


(.MHz) 

(dB) 

(dB) 

MATCHED 

.MISSES 

MISSES 


3.2 

-38.2 

-59.1 

3.0 

0.0 

93.00 

3.1 

4.5 

-30.6 

-49.4 

99.0 

38.0 

150.00 

39.8 

6.8 

-1.9 

-27.1 

329.0 

42.0 

70.00 

82.5 

6.9 

■Bi 

-17.5 

175.0 

20.0 

116.00 

60.1 

9.9 

11.4 

-10.7 

157.0 

54.0 

152.00 

50.8 

10.2 

14.8 

-6.8 

179.0 

73.0 

125.00 

58.9 

13.9 

18.4 

-7.0 

200.0 

47.0 

82.00 

70.9 

14.4 

20.0 

-7.0 

215.0 

58.0 

92.00 

70.0 

17.5 

16.4 

-16.4 

228.0 

66.0 

50.00 

82.0 

IBH 

MM 

-12.6 

254.0 

63.0 

51.00 

83.3 

20.3 

34.0 

m 

188.0 

31.0 

74.00 

71.8 

wmm 

37.6 

6.0 

192.0 

19.0 

68.00 

73.8 

23.2 

43.2 

8.8 

169.0 

20.0 

74.00 

69.5 


PERCENTAGE MATCHED 81.81% 


AVERAGE ERROR -11.10 dB 





























































































Table 16. SUMMER SITE O - 90% ES 


FREQ 

SNR 

ERR 

NUMBER 

ACTUAL 

AMBCO.M 

‘)''o 

(MHz) 

(dB) 

(dB) 

MATCHED 

.MISSES 

.MISSES 


3.2 

-66.3 

-86.8 


0.0 

184.00 

10.2 


-24.7 

-46.2 

41.0 

0.0 

169.00 

19.5 

5.2 

-15.9 

-38.3 

43.0 

0.0 

166.00 

20.6 

6.8 

2.2 

-19.4 

22.0 

0.0 

148.00 

12.9 

6.9 

0.9 

-22.4 

45.0 

0.0 

202.00 

18.2 

9.9 

14.1 

-11.5 

120.0 

0.0 

179.00 

40.1 

10.2 

19.2 

-8.4 

251.0 

0.0 

121.00 

67.5 

13.9 

23.1 

-6.3 

255.0 

4.0 

278.00 

47.8 

14.4 

27.8 

-3.1 

334.0 

7.0 

162.00 

67.3 

17.5 

22.6 

-10.0 

327.0 

17.0 

171.00 

65.7 

18.2 

27.5 

-3.9 

157.0 

17.0 

197.00 

44.4 

2U.3 

33.6 

6.4 

153.0 

27.0 

128.00 

54.4 

20.9 

33.4 

5.3 

81.0 

26.0 

168.00 

32.5 

PERCENTAGE .MATCHED 94.97% AVERAGE ERROR -8.33 dB 
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Table 17. SUMMER SITE O - 50% ES 


FREQ 

SNR 

ERR 

NUMBER 

ACTUAL 

AM BCDM 

“o 

(MUz) 

(dB) 

(dB) 

MATCHED 

MISSES 

MISSES 


3.2 

-66.5 

-87.0 

21.0 

0.0 

184.00 

10.2 

4.9 

-25.2 

-46.7 

41.0 

0.0 

170.00 

19.4 

5.2 

-16.3 

-38.7 

43.0 

0.0 

162.00 

21.0 

6.8 

0.3 

BHI 

22.0 

0.0 

148.00 

12.9 

6.9 

-0.3 

-23.7 

45.0 

0.0 

200.00 

18.4 

9.9 

12.6 

-13.1 

120.0 

0.0 

173.00 

41.0 

10.2 

17.1 

-10.5 

251.0 

0.0 

111.00 

69.3 

13.9 

18.2 

-11.3 

243.0 

16.0 

196.00 

55.4 

14.4 

22.6 

-8.5 

304.0 

37.0 

100.00 

75.2 

17.5 

21.5 

-12.0 

258.0 

86.0 

60.00 

81.1 

18.2 

29.1 

-2.8 

125.0 

49.0 

73.00 

63.1 

20.3 

37.8 

9.3 

124.0 

56.0 


58.8 

20.9 

39.5 

8.9 

55.0 

52.0 

120.00 

31.4 

PERCENTAGE MATCHED 84.80% AVER.AGE ERROR -11.05 dB 
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Table 18. SUMMER SITE O - 0% ES 


FREQ 

SNR 

ERR 

NUMBER 

ACTUAL 

AM BCDM 

O/o 

(MHz) 

(dB) 

(clB) 

.MATCHED 

.MISSES 

.MISSES 


3.2 

-78.9 

-90.4 

21.0 

0.0 

176.00 

10.7 

4.9 

-38.8 

-60.3 

41.0 

0.0 

163.00 

20.1 

5.2 

-28.0 

-50.4 

41.0 

2.0 

132.00 

23.7 

6.8 

-8.4 

-30.0 

21.0 

1,0 

120.00 

14.9 

6.9 

-12.7 

-35.9 

41.0 

4.0 

179.00 

18.6 

9.9 


-13.7 

99.0 

21.0 

135.00 

42.3 

10.2 

14.0 

-13.6 

204.0 

47.0 

97.00 

67.8 

13.9 

20.7 

-8.9 

219.0 

40,0 

221.00 

49.8 

14.4 

24.6 

-6.8 

288.0 

53.0 

108.00 

72.7 

17.5 

21.4 

-12.1 

255.0 

89.0 

64.00 

79,9 

18.2 

30.6 

-0.7 

124.0 

50.0 

66.00 

65.3 

20.3 

37.3 

8.9 

131.0 

49.0 

84.00 

60.9 

20.9 

39.7 

9.0 

56.0 

51.0 

108.00 

34.1 

PERCENTAGE MATCHED 79.11% AVERAGE ERROR -11.87 dB 
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Table 19. SUMMER SITE C - 90% ES 


FREQ 

SNR 

ERR 

NUMBER 

ACTUAL 

AMBC.OM 

°/o 

(MHz) 

(dB) 

(dB) 

MATCHED 

MISSES 

MISSES 


3.2 

-72.2 

-97.5 

15.0 

0.0 

177.00 

7.8 

4.9 

5.5 

-18.4 

12.0 

0.0 

18.00 

40.0 

5.2 

10.6 

-12.8 

23.0 

1.0 

14.00 

62.2 

6.8 

16.9 

1 

o 

52.0 

2.0 

32.00 

61.9 

6.9 

10.7 

-9.8 

42.0 

1.0 

138.00 

23.3 

9.9 

15.7 

-10.4 

59.0 

2.0 

102.00 

36.6 

10.2 

19.7 

-8.7 

215.0 

5.0 

30.00 

87.8 

13.9 

21.8 

-7.8 

192.0 

3.0 

156.00 

55.2 

14.4 

26.1 

-4.6 

230.0 

2.0 

117.00 

66.3 

17.5 

32.3 

-0.6 

139.0 

10.0 

190.00 

42.2 

18.2 

38.4 

6.0 

81.0 

12.0 

133.00 

37.9 

20.3 

44.0 

14.8 

34.0 

5.0 

89.00 

27.6 

20.9 

48.7 

23.9 

2.0 

5.0 

43.00 

on 

PERCENTAGE .MATCHED 95.80% AVERAGE ERROR -6.23 dB 


86 




































































































Table 20. SUMMER SITE C - 50% ES 


FREQ 

SNR 

ERR 

NUMBER 

ACTUAL 

A.MBCOM 


(MHz) 

(dB) 

(dB) 

.MATCHED 

.MISSES 

.MISSES 


3.2 

-72.8 

-98.1 

15.0 

0.0 

170.00 

8.1 

4.9 

4.2 

-19.2 

11.0 

1.0 

27.00 

28.9 

5.2 

9.5 

-13.8 

22.0 

2.0 

25.00 

46.8 

6.8 

15.9 

-8.1 

51.0 

3.0 

38.00 

57.3 

6.9 

9.2 

-11.2 

41.0 

2.0 

145.00 

22.0 

9.9 

14.7 

-11.4 

59.0 

2.0 

106.00 

35.8 

10.2 

17.5 

-10.8 

212.0 

8.0 

48.00 

81.5 

13.9 

17.6 

-12.0 

190.0 

5.0 

119.00 

61.5 

14.4 

22.3 

-8.5 

216.0 

16.0 

67.00 

76.3 

17.5 

31.2 

-2.2 

116.0 

33.0 

70.00 

62.4 

18.2 

39.8 

7.3 

75.0 

18.0 

89.00 

45.7 

20.3 

45.6 

13.5 

29.0 

10.0 

83.00 

25.9 

20.9 

48.7 

23.9 

2.0 

5.0 

42.00 

4.5 

PERCENTAGE .MATCHED 90.82% AVERAGE ERROR -8.88 dB 
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Table 21. SUMMER SITE C - 0% ES 


FREQ 

SNR 

ERR 

NU.MBER 

ACTUAL 

AMBCO.M 


(MHz) 

(dB) 

(dB) 

.MATCHED 

MISSES 

MISSES 


3.2 

-83.7 

-109.1 

15.0 

0.0 

169.00 

8.2 

Bam 

-12.2 

-35.6 

11.0 

1.0 

47.00 

19.0 

5.2 

-3.5 

-26.7 

22.0 

2.0 

24.00 

47.8 

6.8 

4.4 

-19.5 

48.0 

6.0 

40.00 

54.5 

6.9 

-1.3 

-21.8 

41.0 

2.0 

135.00 

23.3 

9.9 

5.8 

-20.1 

50.0 

11.0 

115.00 

30.3 

10.2 

9.1 

-19.3 

178.0 

42.0 

72.00 

71.2 

13.9 

14.67 

-15.0 

166.0 

29.0 

126.00 

56.8 

14.4 

19.6 

-11.3 

196.0 

36.0 

95.00 

67.4 

17.5 

32.9 

-0.7 

105.0 

44.0 

82.00 

56.1 

18.2 

42.3 

9.9 

73.0 

20.0 

78.00 

48.3 

20.3 

46.9 

14.8 

29.0 

10.0 

68.00 

29.9 

20.9 

48.7 

23.9 

2.0 

5.0 

42.00 

4.5 

PERCENTAGE .MATCHED 81.82% AVERAGE ERROR -13.33 dB 
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Table 22. SUMMER SITE D - 90% ES 


FREQ 

SNR 

ERR 

NUMBER 

ACTUAL 

AMBCO.M 

O/o 

(MHz) 

(dB) 

(dB) 

MATCHED 

MISSES 

MISSES 


3.2 

-4.2 

-36.3 

36.0 

7.0 

22.00 

62.1 

4.9 

5.1 

-19.6 

109.0 

0.0 

44.00 

71.2 

5.2 

10.8 

-18.3 

128.0 

0.0 

36.00 

78.0 

6.8 

16.2 

-10.0 

178.0 

0.0 

69.00 

72.1 

6.9 

21.9 

-5.3 

133.0 

0.0 

94.00 

58.6 

9.9 

15.6 

-12.4 

164.0 

0.0 

175.00 

48.4 

10.2 

20.5 

-11.1 

359.0 

0.0 

90.00 

80.0 

13.9 

28.4 

-8.7 

198.0 

1.0 

206.00 

49.0 

14.4 

32.8 

-5.1 

340.0 

0.0 

138.00 

71.1 

17.5 

42.6 

1.6 

231.0 

44.0 

196.00 

54.1 

18.2 

48.2 

9.3 

74.0 

55.0 

189.00 

28.1 

20.3 

48.2 

10.0 

46.0 

107.0 

98.00 

31.9 

20.9 

47.0 

11.4 

15.0 

31.0 

125.00 

10.7 

PERCENTAGE MATCHED 89.14% AVERAGE ERROR 

-8.02 dB 
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Table 23. SUMMER SITE D - 50% ES 


FREQ 

SNR 

ERR 

NUMBER 

ACTUAL 

AMBCOM 

«/o 

(MHz) 

(dB) 

(dB) 

MATCHED 

.MISSES 

MISSES 


3.2 

-4.6 

-36.7 

36.0 

7.0 

22.00 

62.1 

4.9 

4.4 

-20.4 

109.0 

0.0 

46.00 

70.3 

5.2 

10.1 

-18.9 

128.0 

0.0 

35.00 

78.5 

6.8 

15.2 

-11.1 

178.0 

0.0 

67.00 

72.7 

6.9 

21.0 

-6.3 

133.0 

0.0 

101.00 

56.8 

9.9 

14.1 

-14.0 

163.0 

1.0 

167.00 

49.4 

10.2 

17.6 

-14.1 

356.0 

3.0 

63.00 

85.0 

13.9 

25.1 

-13.0 

171.0 

28.0 

103.00 

62.4 

14.4 

27.6 

-11.4 

273.0 

67.0 

64.00 

81.0 

17.5 

43.1 

0.2 

184.0 

91.0 

64.00 

74.2 

18.2 

49.0 

8.8 

53.0 

76.0 

53.00 

:o.o 

20.3 

49.8 

5.3 

26.0 

127.0 

8.00 

76.5 

20.9 

46.9 

3.8 

5.0 

41.0 

OO 

38.5 


PERCENTAGE MATCHED 80.45% AVERAGE ERROR -11.44 dB 
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Table 24. SUMMER SITE D - 0% ES 


FREQ 

SNR 

ERR 

NUMBER 

ACTUAL 

AMBCOM 

°/o 

(.MHz) 

(dB) 

(dB) 

MATCHED 

MISSES 

MISSES 


3.2 

-14.7 

-46.4 

34.0 

9.0 

25.00 

57.6 

4.9 

-13.1 

-37.8 

100.0 

9.0 

84.00 

54.3 

5.2 

-5.2 

-34.5 

121.0 

7.0 

60.00 

66.9 

6.8 

8.0 

-18.5 

163.0 

15.0 

40.00 

80.3 

6.9 

12.1 

-15.5 

115.0 

18.0 

69.00 

62.5 

9.9 

9.5 

-18.9 

144.0 

20.0 

96.00 

60.0 

10.2 

12.9 

-19.1 

309.0 

50.0 

56.00 

84.7 

13.9 

20.4 

-17.5 

173.0 

26.0 

109.00 

61.3 

14.4 

24.9 

-13.7 

279.0 

61.0 

59.00 

82.5 

17.5 

43.4 

0.5 

182.0 

93.0 

64.00 

74.0 

18.2 

48.6 

8.6 

60.0 

69.0 

60.00 

50.0 

20.3 

47.5 

4.8 

29.0 

124.0 


70.7 

20.9 

46.9 

3.8 

5.0 

41.0 

10.00 

33.3 

PERCENTAGE .MATCHED 75.98% AVERAGE ERROR -16.96 dB 


\ I. 
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APPENDIX D. COMPARIABLE STATISTICS 


The tables and figures listed in this Appendix were compiled from a subset of 
AMBCOM data. Earlier computer-based model validations [Ref. 14 and 13] modeled 
six of the 13 Project NONCENTRIC frequencies, e.g., 6.8, 6.9, 9.9, 13.9, 17.5 and 20.3 
.MHz. In addition, only site O data were examined. The data presented in this Appendix 
are for those frequencies and site O only. The figures are directly compareable to the 
figures produced in the earlier studies. 


Table 25. WINTER SITE O STATISTICS (NORMALIZED TO ZERO MEAN 
ERROR). 


% Sporadic E 

90% 

50‘>/o 

0% 

Standard Dev. (dB) 

15.89 

17.87 

19.86 

Median Error (dB) 

-0.1 

0.9 

5.8 

# of Samples 

1012 

928 

750 

% .Matched 

89.80 

82.34 

66.55 


Table 26. SUMMER SITE O STATISTICS (NORMALIZED TO ZERO MEAN 


ERROR). 


% Sporadic E 

90% 

50°/o 

O'^.o 

Standard Dev. (dB) 

14.16 

15.38 

17.52 

Median Error (dB) 

0.2 

-0.6 

1.1 

^ of Samples 

922 

812 

766 

"/o Matched 

95.05 

83.71 

78.97 
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Figure 33. The 1989 winter campaign average error for all three sporadic E models. 
(lU dB bins) 
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Figure 34. The 1989 winter campaign cumulative error of absolute value of the av¬ 
erage error. 
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Figure 35. The 1988 suiniiier campaign average error for all tliree sporadic E mod¬ 
els. (10 dB bins) 
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Figure 36. The 1989 summer campaign cumulative error of absolute value of the 
average error. 


96 


ABSOLUTE- VALUES OF ERROR (dB) 







LIST OF REFERENCES 


1. J.K. Hargreaves, "The high latitude ionosphere: dynamical aspects and models", 
paper presented at the NATO AGARD Lecture Series, Media EfTects on Electronic 
Systems in the High Latitude Region, LS-162, September I9S8. 

2. Interview with Tudor Jones, Professor, University of Leicester, July 1990. 

3. W. D. Parkinson, Introduction to Geomagnetism, Scottish Academic Press Ltd., 
1983. 

4. Kenneth Davies, Ionospheric Radio,?cXQX Peregrinus Ltd., London, United 
Kingdom, 1989. 

5. Helen E. Coffey, "Geomagnetic and Solar Data", Journal of Geophysical Research, 
Vol. 94, No A6, p. 7020, June 1989. 

6. J.K. Hargreaves, "The high-latitude ionosphere: geophysical basis", paper presented 
at the NATO AGARD Lecture Series, Media Effects on Electronic Systems in the 
High Latitude Region, LS-162, September 1988. 

7. R. D. Hunsucker and H. F. Bales, "Survey of polar and auroral effects on HF 
propagation", Radio Science, Vol. 4, No. 4, pp. 347-365, April 1969. 

8. R. R. Vondrak and K. L. Miller, "A high-latitude phenomenological model of 
auroral precipitation and ionospheric effects". Radio Science, Vol. 20, No. 3, pp. 
431-438, May-June 1985. 

9. R. R. Vondrak and others, "Chatanika model of the high-latitude ionosphere for 
application to HF propagation prediction", RADC-TR-78-7, Menlo Park, CA, 
January 1978. 


97 






10. J. K. Hargreaves, "Auroral absorption of HF radio waves in the ionosphere: a re¬ 
view of results from the first decade of riometr>’", Proceedings of the IEEE, Vol. 57, 

No.8, August 1969. 

11. D.P. Roesler, Analysis of Data from A Transauroral HE Experiment, Report 
WP88-2010, Rockwell International, March 1988. 

12. Nikhil Dave', "Comparison with oblique sounder data of high- latitude HF pred¬ 
ictions from RADAR C and A.MBCOM computer programs", paper presented at 
the Ionospheric Effects Symposium, Springfield, VA, May 1987. 

13. Marcos D. Tsolekas, "A comparison of high-latitude ionospheric propagation 
predictions from lONCAP-PC 2.5 with measured data," Master's thesis, Naval 
Postgraduate School, Monterey, CA, December 1990. 

14. Stefanos S. Gikas, "A comparison of high-latitude ionospheric propagation pred¬ 
ictions from Advanced PROPHET 4.0 with measured data," Master's thesis. Naval • 
Postgraduate School, Monterey, CA, December 1990. 

15. V. Elaine Hatfield and Georgellen Smith, AMBCOM USER'S GUIDE FOR EN¬ 
GINEERS, SRI International, 1987. 

16. V. Elaine Hatfield and others, AMBCOM USER'S GUIDE FOR PROGRAM¬ 
MERS: Part One — Description, SRI International, 1987. 

17. V. Elaine Hatfield and others, AMBCOM USER'S GUIDE FOR PROGRAM¬ 
MERS: Part Two — Listings, SRI International, 1987. 

18. John Griffiths, Radio Wave Propagation and Antennas, Prentiss-Hall International, 
Englewood Cliffs, NJ, 1987. 

19. Helen E. Coffey, "Geomagnetic and Solar Data", Journal of Geophysical Research, 

Vol.94, No. A7, p. 9161, July 1989. 


98 








20. Helen E.CofTey, "Geomagnetic and Solar Data", Journal of Geophysical Research, 
Vol.93, No. A12, p. 14732, December 1988. 

21. Helen E. Colley, "Geomagnetic and Solar Data", Journal of Geophysical Research, 
Vol.94. No. Al, p. 511, 1 January 1989. 

22. T. R. Hartz and N. M. Brice,"The general pattern of auroral particle precipitation," 
Planetary and Space Science, Vol. 15.1, p. 301-327, 1967. 

23. Sonceniric:Automatic Signal Recognition, Ionospheric Physics Group, Leicester 
University, December 1989. 

24. Telephone conversation between V. E. Hatfield, SRI, International and the author, 
25 February' 1991. 


99 






